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ABSTRACT 

The  successive  stages  of  the  fireball  due  to  a  nuclear  explosion  in 

air  are  defined  (Sec.  2).  This  paper  is  chiefly  concerned  with  Stage  C, 

from  the  minimum  in  the  apparent  fireball  temperature  to  the  point  where 

the  fireball  becomes  transparent.  In  the  first  part  of  this  stage  (C  I), 

the  shock  (which  previously  was  opaque)  becomes  transparent  due  to 

decreasing  pressure.  The  radiation  comes  from  a  region  in  which  the 

temperature  distribution  is  given  essentially  by  the  Taylor  solution;  the 

radiating  layer  is  given  by  the  condition  that  the  mean  free  path  is 

about  l/50  of  the  radius  (Sec.  3).  The  radiating  temperature  during  this 

-l/b 

stage  increases  about  as  p  '  ,  where  p  is  the  pressure. 

To  supply  the  energy  for  the  radiation,  a  cooling  wave  proceeds  from 
the  outside  into  the  hot  interior  (Sec.  5) •  When  this  wave  reaches  the 
isothermal  sphere,  the  temperature  is  close  to  its  second  maximum.  There¬ 
after,  the  character  of  the  solution  changes;  it  is  now  dominated  by  the 
cooling  wave  (Stage  C  II).  The  temperature  would  decrease  slowly  (as 
p"^)  if  the  problem  were  one -dimensional,  but  in  fact  it  is  probably 
nearly  constant  for  the  three -dimensional  case  (Sec.  6).  The  radiating 
surface  shrinks  slowly.  The  cooling  wave  eats  into  the  isothermal  sphere 
until  this  is  completely  used  up.  The  inner  part  of  the  isothermal  sphere, 
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i.e.,  the  part  which  has  not  yet  been  reached  by  the  cooling  wave,  con¬ 
tinues  to  expand  adiab at i c ally;  it  therefore  cools  very  slowly  and 
remains  opaque. 

After  the  entire  isothermal  sphere  is  used  up,  the  fireball  becomes 
transparent  and  the  radiation  drops  rapidly.  The  ball  will  therefore  be 
left  at  a  rather  high  temperature  (Sec.  7)>  about  5000°. 

The  cooling  wave  reaches  the  isothermal  sphere  at  a  definite  pres¬ 
sure  pc  «  5(Pl/p0)l/5  bars,  where  is  the  ambient  and  pQ  the  sea  level 
density.  The  radiating  temperature  at  this  time  is  about  10,000°.  The 
slight  dependence  of  physical  properties  on  yield  is  exhibited  in  approx¬ 
imate  formulae. 
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1.  INTRODUCTION 


Tne  radiation  from  the  fireball  has  been  studied  intensively  by 
many  authors.  Already  in  the  Simmer  Study  at  Berkeley  in  1942,  Bethe 
and  Teller"*"  found  that  the  energy  transmitted  by  a  nuclear  explosion 
into  air  is  immediately  converted  into  x-rays,  and  studied  the  qualita¬ 
tive  features  of  the  transmission  of  these  x-rays.  At  Los  Alamos, 

2 

Marsnak  and  others  showed  that  this  radiation  propagates  as  a  wave, 

3 

with  a  sharp  front.  Hirschf elder  and  Magee  gave  the  first  comprehensive 
treatment  of  this  early  phase  of  the  fireball  development,  and  also 
studied  some  of  the  later  phases,  especially  the  role  of  NOg. 

Many  optical  observations  have  been  made  in  the  numerous  tests  of 
atomic  weapons.  Some  of  the  results  are  contained  in  "Effects  of  Nuclear 
Weapons, pp.  70-84  (see  also  pp.  316-368).  A  summary  of  the  spectro¬ 
scopic  observations  up  to  1956  was  compiled  by  DeWitt.^  Careful  scrutiny 
of  the  extensive  observational  material  would  undoubtedly  give  a  wealth 
of  further  information. 

On  the  theoretical  side,  there  has  been  some  analytical  and  a  good 
deal  of  numerical  work.  Analytical  work  has  concentrated  on  the  early 
phases.  One  of  the  most  recent  analytical  papers  on  the  early  flow  of 
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radiation  (Stage  A  of  Sec.  2)  is  by  Freeman.  Brode  and  Gilmore  treat 
also  Stage  B,  the  radiation  from  the  shock  front,  with  particular  emphasis 
on  the  dependence  on  altitude. 

O 

The  most  complete  numerical  calculation  has  been  done  by  Brode  on 
a  sea  level  megaton  explosion.  We  shall  use  his  results  extensively, 
but  for  convenience  we  shall  translate  them  to  a  yield  of  1  megaton. 
Wherever  the  phenomena  are  purely  hydrodynamic,  we  may  simply  scale  the 
linear  dimensions  and  the  time  by  the  cube  root  of  the  yield,  and  this 
is  the  principal  use  we  shall  make  of  Brode's  results,  e.g.  in  Sec.  5b. 
Where  radiation  is  important,  this  scaling  will  give  only  a  rough  guide. 
Brode  calculates  pressures,  temperatures,  densities,  etc.,  as  functions 
of  time  and  radius,  for  scaled  ( 1-megaton)  times  of  about  10 to  10 
seconds.  The  calculations  show  clearly  the  stages  in  fireball  and  shock 
development ,  as  defined  in  Sec.  2,  at  least  Stages  A  to  C. 

Q 

Brode  and  Meyerott  have  considered  the  physical  phenomena  involved 
in  the  optical  "opening"  of  the  shock  after  the  minimum  of  radiation. 

Stage  C  I  in  the  nomenclature  of  Sec.  2,  especially  the  decrease  in 
opacity  due  to  decreased  density  and  to  the  dissociation  of 

Zel’dovich,  Kompaneets  and  Raizer"*"^  have  investigated  how  the  energy 
for  the  radiation  is  supplied  after  the  radiation  minimum  and  have  intro¬ 
duced  the  concept  of  a  "cooling  wave"  moving  into  the  hot  fireball.  The 
present  report  is  largely  concerned  vith  an  extension  of  the  ideas  of 
Zel'dovich  et  al.  to  the  actual  case  of  density  varying  vith  time,  more 
general  opacity  function,  radiation  absorption  varying  vith  vave  length, 
etc. 
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Much  effort  has  been  devoted  to  the  calculation  of  fireballs  at 

various  altitudes.  Brode"1"1  made  such  calculations  in  1958,  in  connection 

12 

vith  the  test  series  of  that  year.  Gilmore  made  a  prediction  of  the 
Bluegill  explosion  in  19^2 .  Since  then,  many  more  refined  calculations 
of  Bluegill  have  been  made. 

For  any  understanding  of  fireball  phenomena  it  is  essential  to  have 

a  good  equation  of  state  for  air  and  good  tables  of  absorption  coeffi- 

13 

cients.  For  the  equation  of  state,  we  have  used  Gilmore's  tables, 

14 

although  Hilsenrath *s  give  more  detail  in  some  respects.  The  two  cal¬ 
culations  agree.  Gilmore's  equation  was  approximated  analytically  by 

„  „  8 

Brode . 

For  the  absorption  coefficient  we  have  used  the  tables  by  Meyerott 

et  al.,1^1^  which  extend  to  12,000°.  These  were  supplemented  by  the 

17  18 
work  of  Kivel  and  Bailey  and  more  recent  work  by  Taylor  and  Kivel  on 

the  free-free  electron  transitions  in  the  field  of  neutral  atoms  and 

molecules.  At  higher  temperatures,  there  are  calculations  by  Gilmore 

19  20  21 

and  Latter,  ^  Karzas  and  Latter,  and  curves  by  Gilmore  which  are 

brought  up  to  date  periodically.  The  most  recent  calculation  on  the 

absorption  of  air  at  about  18,000°  and  above  have  been  done  by  Stuart 

op 

and  Pyatt.  This  temperature  range  is  not  of  great  importance  for  the 
problem  of  this  paper,  but  is  important  for  the  expansion  of  the  iso¬ 
thermal  sphere  inside  the  shock  wave  before  it  is  reached  by  the  cooling 
wave. 

Q 

Brode  has  used  the  average  of  the  absorption  coefficient  over 
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frequency,  the  opacity,  which  is  sufficient  for  treating  the  internal 
flow  of  radiation.  A  realistic  treatment  of  the  flow  to  the  outside 
requires  the  absorption  coefficient  as  a  function  of  frequency;  Brode 
merely  wanted  to  obtain  reasonable  overall  results  for  this  flow.  He 
approximated  the  opacity  by  an  analytic  expression.  Also  in  most  of  the 
other  work  cited  above  an  average  opacity  has  been  used.  An  exception 
is  some  of  the  recent  work  on  the  radiation  flow  in  high  altitude  explo¬ 
sions  (Bluegill),  where  the  frequency  dependence  must  be,  and  has  been, 

21 

taken  into  account .  Gilmore  has  calculated  and  made  available  curves 
of  effective  opacity,  in  which  the  radiation  mean  free  path  was  averaged 
(using  a  Rosseland  weighting  factor)  only  over  those  frequencies  for 
which  it  is  less  than  1  kilometer. 

This  list  of  references  on  work  on  the  fireball  dynamics  and  opac¬ 
ity  is  far  from  complete. 


2.  PHASES  IN  DEVELOPMENT 

The  energy  from  a  nuclear  explosion  is  transmitted  through  the  outer 
parts  of  the  weapon,  including  its  case,  either  by  radiation  (x-rays)  or 
by  shock  or  both.  Whichever  the  mode  of  transmission  inside  the  weapon, 
once  the  energy  gets  into  the  surrounding  air,  the  energy  will  be  trans¬ 
ported  by  x-rays.  This  is  because  the  air  will  be  heated  to  such  a  high 
temperature  (a  million  degrees  or  more)  that  transport  by  x-rays  is  much 
faster  than  by  hydrodynamics.  This  stage  of  energy  transport  (Stage  A) 
has  been  extensively  studied  by  many  authors  (e.g.,  Hirschfelder  and 
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Magee  in  Report  LA -2000)  and  will  therefore  not  he  further  considered 
here. 

During  Stage  A,  temperatures  are  very  high.  The  Planck  spectrum  of 
the  air  is  in  the  x-ray  or  far  ultraviolet  region,  and  hence  is  immedi¬ 
ately  absorbed  by  the  surrounding  air.  The  very  hot  air  is  therefore 
surrounded  by  a  cooler  envelope,  and  only  this  envelope  is  visible  to 
observers  at  a  distance.  The  observable  temperature  therefore  has  little 
physical  significance.  It  is  observed  that  the  size  of  the  luminous 
sphere  increases  rapidly,  and  the  total  emission  also  increases,  up  to 
a  first  maximum. 

When  the  temperature  of  the  central  sphere  of  air  has  fallen,  by 
successive  emission  and  re-absorption  of  x-rays,  to  about  300,000°,  a 
hydrodynamic  shock  forms.  The  shock  now  moves  faster  than  the  tempera¬ 
ture  could  propagate  by  radiation  transport.  The  shock  therefore  sepa¬ 
rates  from  the  very  hot,  nearly  isothermal  sphere  at  the  center.  This 
is  Stage  B  in  the  development.  The  shock  moves  by  simple  hydrodynamics. 
Its  front  obeys  the  Hugoniot  relations,  the  density  being  given  by  (3«*0 • 
Behind  the  front,  the  air  expands  adiabatically,  and  at  a  radius  of  80# 
of  the  shock  radius,  the  density  is  apt  to  have  fallen  by  a  factor  of  10 
or  more  compared  to  the  shock  density,  while  the  temperature  has  risen 
by  a  comparable  factor,  (3*l6).  Thus  the  interior  is  at  very  low  density, 
and  hence  the  pressure  must  be  nearly  constant  (otherwise  there  would  be 
very  large  accelerations  which  soon  would  equalize  the  pressure).  This 
greatly  simplifies  the  structure  of  the  shock,  and  leads  to  such  simple 
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relations  as  (5*2)  between  shock  radius  and  pressure. 

Well  inside  the  shock,  the  "isothermal  sphere"  pursues  its  separate 
history.  It  continues  to  engulf  more  material  because  radiation  flow 
continues,  though  at  a  reduced  rate.  H.  Brode  has  kindly  calculated  for 
me  the  temperature  histories  of  several  material  points,  based  on  his 
paper  PM -2248.  These  histories  clearly  exhibit  the  expansion  of  the  iso¬ 
thermal  sphere  in  material  coordinates.  The  expansion  can  also  be  treated 
by  a  semi -analytic  method  which  I  hope  to  discuss  in  a  subsequent  paper. 
The  isothermal  sphere  remains  isolated  from  the  outside  world  until  it 
is  reached  by  the  cooling  wave.  Sec.  5* 

The  radiation  to  the  outside  now  comes  from  the  shock.  Early  in 
Stage  B,  the  shock  has  a  precursor  of  lower  temperature,  caused  by  ultra¬ 
violet  radiation  from  the  shock,  and  the  observable  temperature  is  still 
below  the  shock  temperature  (Stage  B  I).  However,  the  observable  radius 
is  very  near  the  shock  radius.  Later,  as  the  shock  front  cools  down, 
the  shock  radiates  directly  and  its  temperature  becomes  directly  observ¬ 
able  (Stage  B  II )  •  The  first  maximum  in  visible  radiation  probably 
occurs  between  Stages  B  I  and  B  II.  As  the  shock  cools  down,  the  radi¬ 
ation  from  the  shock  front  decreases,  and  the  observable  temperature 
decreases  to  a  minimum  (Ref.  4,  par.  2.113,  p.  75)  of  about  2000°. 

When  the  shock  is  sufficiently  cool,  its  front  becomes  transparent, 
and  one  can  look  into  it  to  higher  temperatures  (Stage  C).  The  central 
isothermal  sphere,  however,  remains  opaque  and,  for  some  time,  invisible. 
Because  higher  temperatures  are  now  revealed,  the  total  radiation  increases 
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toward  a  second  maximum.  This  stage  has  been  very  little  considered 
theoretically,  except  in  numerical  calculations,  and  forms  the  subject 
of  this  report. 

We  shall,,  show  that  for  some  time  in  Stage  C,  the  radiation  comes 
from  the  air  between  isothermal  sphere  and  shock  front  (Stage  C  I). 
During  this  time,  the  radiation  can  be  calculated  essentially  from  the 
temperature  distribution  which  is  set  up  by  the  adiabatic  expansion  of 
the  material  behind  the  shock  (Secs.  3  and  5f  ,g) •  The  temperature  and 
total  intensity  of  the  radiation  increase  with  time  toward  the  larger, 
second  maximum. 

The  energy  for  the  radiation  is  largely  supplied  by  a  cooling  wave 
(Sec.  5)  which  gradually  eats  into  the  hot  interior.  When  this  cooling 
wave  reaches  the  isothermal  sphere,  the  radiation  temperature  reaches 
its  maximum  (Sec.  5e);  it  then  declines  again  as  the  cooling  wave  eats 
more  deeply  into  the  isothermal  sphere  (Stage  C  II).  This  process  ends 
by  the  isothermal  sphere  being  completely  eaten  away  (Sec.  6). 

After  this  has  happened,  the  entire  fireball.,  isothermal  sphere  and 
cooler  envelope,  is  transparent  to  its  own  thermal  radiation  (Stage  D). 
The  molecular  bands,  which  previously  appeared  in  absorption,  now  appear 
in  emission  (Stage  D  I).  Emission  will  lead  to  further  cooling  of  the 
fireball,  though  more  slowly  than  before.  Soon,  when  the  temperature 
falls  below  about  6000°,  the  emission  becomes  very  weak,  and  subsequent 
cooling  is  almost  entirely  adiabatic  (Stage  D  IX ).  At  sea  level,  the 
pressure  may  go  down  to  1  bar  before  the  temperature  falls  below  6000° ; 
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in  this  case  there  is  no  Stage  D  II.  At  higher  elevation,  there  usually 
is. 

The  fireball  will  then  remain  hot,  at  about  6000°  or  a  somewhat 
lower  temperature  due  to  adiabatic  expansion  in  Stage  D  II.  The  only 
process  which  can  now  lead  to  further  cooling  is  the  rise  of  the  fire¬ 
ball,  which  leads  to  further  adiabatic  expansion  and,  more  important,  to 
turbulent  mixing  at  the  surface  with  the  ambient  air  (Stage  E).  The 
time  required  for  this  is  typically  10  seconds  or  more,  being  determined 
by  buoyancy. 

At  very  high  altitude,  the  shock  wave  never  plays  an  important  part, 
but  radiation  transport  continues  until  the  temperature  gets  too  low  for 
effective  emission.  In  other  words,  Stage  A  continues  to  the  end.  Of 
course,  a  shock  does  form,  but  it  is,  so  to  speak,  an  afterthought,  and 
it  plays  little  part  in  the  distribution  of  energy.  At  medium  altitude, 
let  us  say,  10  to  30  kilometers,  the  stages  are  much  the  same  as  at  sea 
level  but  the  shock  wave  becomes  transparent  earlier,  i.e.,  at  a  higher 
temperature,  because  the  density  is  lower;  this  means  that  the  minimum 
emission  comes  earlier.  Stage  C  proceeds  similarly  to  sea  level,  but 
at  the  second  maximum  of  radiation  the  pressure  is  still  much  above 
ambient,  therefore  Stage  C  II  involves  a  greater  radial  expansion  of  the 
isothermal  sphere  than  at  sea  level  which  proceeds  simultaneously  with 
the  inward  motion  of  the  cooling  wave.  Moreover,  there  is  much  adia¬ 
batic  expansion  after  the  cooling  wave  has  penetrated  to  the  center. 

The  temperature  at  which  radiation  stops  is  higher,  due  to  the  lower 
density. 
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We  believe  that  the  theory  developed  in  this  paper  will  be  useful 
in  studying  the  dependence  of  phenomena  on  altitude  (ambient  pressure), 
but  we  have  not  yet  exploited  it  for  this  purpose. 

3.  RADIATION  FROM  BEHIND  THE  SHOCK  (Stage  C  I) 

a.  Role  of  NOg 

The  diatomic  species  in  equilibrium  air,  both  neutrals  and  ions, 
show  very  little  absorption  in  the  visible  at  temperatures  up  to  about 
4000°K.  This  is  shown  clearly  in  the  tables  by  Meyerott  et  al.1^  We 
define  "the  visible"  for  the  purposes  of  this  paper,  arbitrarily  and 
incorrectly,  as  the  frequency  range 

hv  =  l/2  to  2-3 A  ev 

=  4050  to  22,300  cm"1 

=  2.48  to  0.45  p  (3.1) 

Then,  even  at  a  density  as  high  as  lOp^  (Pq  =  density  of  air  at  NTP 
-3  3 

=  1.29  x  10  gm/cnr  )>  the  mean  free  path  is  never  less  than  100  meters 
at  4000°,  1000  meters  at  3000°,  and  still  longer  at  lower  temperature. 
These  values  refer  to  hv  =  2-5/8  ev;  for  lower  frequencies,  the  mean 

free  path  is  even  longer. 

In  the  sea-level  shock  wave  from  a  megaton  explosion,  the  temper¬ 
ature  range  from  3000  to  4000°  occupies  a  distance  of  about  10  meters 
(see  Sec.  3b).  Thus  this  region  is  definitely  transparent,  even  at  the 
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highest  possible  density  of  about  lOp^.  For  explosions  at  higher  alti¬ 
tude,  this  conclusion  is  even  more  true. 

The  tables  by  Meyerott  et  al.  do  not  include  absorption  by  triatomic 
(and  more  complicated)  molecules.  Of  these,  N0g  is  known  to  have  strong 
absorption  bands  in  the  visible.  After  this  paper  vas  completed,  I 
received  new  calculations  by  Gilmore  '  which  include  the  effect  of  N0g. 
The  effect  is  very  striking  as  is  shown  by  Table  I,  which  gives  the 
absorption  coefficient  for  the  "typical"  frequency  hv  =  2-l/8  ev,  and 
for  several  temperatures  and  densities  (the  absorption  is  strong  from 

about  1-3/4  to  2-3/4  ev). 


Table  I.  Absorption  Coefficients  at  hi/  =  2-1/8  ev 
with  and  without  N0g 
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3.1-2 
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Wote:  For  each  value,  the  power  of  10  is  indicated  by  a  superscript. 

Because  HOg  is  triatomic,  its  absorption  depends  strongly  on  density. 

As  the  shocked  gas  ejqpands,  the  N0g  dissociates  and  the  gas  becomes  trans- 

g 

parent.  Brode  and  Meyerott  have  calculated,  under  reasonable  assump¬ 
tions,  the  effect  of  this  dissociation  on  the  optical  properties  of  the 
fireball.  We  shall  not  discuss  the  effect  of  N0g  any  further,  but  shall 
assume  that  this  substance  has  almost  disappeared  by  the  time  we  are  con¬ 
sidering. 
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b.  Temperature  Distribution  behind  the  Shock 


We  wish  to  calculate  the  temperature  distribution  behind  the  shock. 
We  can  do  this  because  the  material  which  has  gone  through  the  shock 
expands  very  nearly  adiabatically,  as  long  as  it  is  not  engulfed  by  the 
internal,  hot  isothermal  sphere.  We  are  interested  in  the  period  when 
the  shock  temperature  goes  from  about  4000  to  a  few  hundred  degrees,  i.e., 
until  the  strong  cooling  wave  (Sec.  5)  starts.  For  a  1-megaton  explosion, 
this  corresponds  roughly  to  t  =  0.05  to  0.25  second. 

At  a  given  time,  the  pressure  is  nearly  constant  over  most  of  the 
volume  inside  the  shock,  except  for  the  Immediate  neighborhood  of  the 
shock;  the  shock  pressure  is  roughly  twice  this  constant,  inside  pressure. 
Coinparing  two  material  elements  in  the  "inside"  region,  we  may  calculate 
their  relative  temperatures  if  we  know  their  temperatures  when  the  shock 
traversed  them,  and  assume  adiabatic  expansion  from  there  on. 

A  material  element  which  is  initially  at  point  r  will  be  shocked 
when  the  shock  radius  is  r.  The  shock  pressure  at  this  time  is* 

Pg(r)  =  20 (r^v  -  l)Yr~5  (5-2) 

where 

Y  =  yield  in  megatons 
r  =  radius  in  kilometers 

y'  .  i  =  L.  =  _ pressure _  .  . 

'  pE  energy  per  unit  volume 

dotations  for  thermodynamic  quantities  similar  to  those  of  Gilmore 
(Report  RM-1545). 
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and  the  average  of  y *  is  taken  over  the  volume  inside  the  shock  wave. 
The  basis  of  (3.2)  is  that  the  total  energy  in  the  shocked  volume,  Y, 
is  the  volume  times  the  average  energy  per  -unit  volume,  the  latter  is 
the  average  pressure  divided  by  7'  -  1,  and  the  average  pressure  is 
close  to  one -half  of  the  shock  pressure. 

The  density  at  the  shock  is 


7'  +  1 

s 


(5-^) 


where  the  subscript  5  refers  to  shock  conditions.  Now  an  examination  of 
Gilmore's  tables  shows  that  7'  does  not  change  very  much  along  an  adi- 
abat.  As  an  example,  we  list  in  Table  II  certain  quantities  referring 
to  some  adiabats  which  will  be  particularly  important  for  our  theory. 
These  are  the  adiabats  for  which  the  temperature  T  is  between  4-000  and 
12,000°  at  a  density  of  0.1pQ.  In  the  second  line  we  list  the  temper¬ 
ature  T  for  the  same  entropy  S  at  a  density  p  =  10pA.  This  is  close 

enough  to  the  shock  density  (3»4-)  so  we  may  consider  T  as  the  temper- 

s 

ature  of  the  same  material  when  the  shock  wave  went  through  it.  (Adia¬ 
batic  expansion,  i.e.,  no  radiation  transport,  is  assumed.)  The  third 
line  gives  7*  -  1  for  the  "present"  conditions,  p  =  0.1pQ  and  T,  the 
fourth  line  is  the  same  quantity  for  the  shock  conditions.  It  is  evident 
that  7*  -  1  is  nearly  constant  for  T  =  4-000  and  6000°,  not  so  constant 
for  8000  and  12,000°.  On  the  average  7'  •  1  «  0.l8.  The  last  two  lines 
in  Table  II  give  the  number  of  particles  (atoms,  ions,  electrons)  per 
original  air  molecule  under  "present"  and  shock  conditions. 
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Table  II.  Adiabats 
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1.6 

2.06 

3 

Assuming 

an  adiabat  of 

constant  7 

=7’,  the  density  of  a  mass 

element  is 


(3-5) 


Now  p  is  a  constant,  and  at  any  given  time,  p  is  the  same  for  al 1  mass 
s 

elements  except  those  very  close  to  the  shock,  hence 


p  ~  p 

K  ^s 

(~  means  proportional  to) . 

If  we  now  introduce  the  abbreviation 


(3.6) 


(3-7) 


which  is  proportional  to  the  mass  inside  the  mass  element  considered, 
and  if  we  use  (3*2),  we  find  that  at  any  given  time 
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(3.8) 


p  ~  m 


The  radius  R  is  given  by 


■JH 


dm  7-1  , 

-577  =  m  - const. 

m  ' '  1 


(3.9) 


We  shall  set  the  constant  equal  to  zero  which  amounts  to  the  (incorrect) 
assumption  that  (3*8)  holds  down  to  m  =  0.  Actually,  the  isothermal 
sphere  gives  the  constant  a  finite,  positive  value. 

To  find  the  temperature  distribution,  we  note  that  the  enthalpy 


7-1  P 


(3.10) 


We  are  using  the  enthalpy  rather  than  the  internal  energy  because  the 
interior  of  the  shock  is  at  constant  pressure,  not  constant  density. 
The  thermodynamic  equation  for  H  is 


T  dS  =  dH  -  v  dp 


(3.11) 


At  given  pressure,  i.e.,  given  time,  (3*8)  to  (3*10)  give 


H  ~  p"1  ~  m"1//7  ~  R"3/(7-l) 


(3.12) 
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New  the  approximate  relation  between  internal  energy  and  temperature  is 
13 

given  by  Gilmore,  Fig.  5,  viz.,* 

E  =  4.25  x  10n  T*1^  erg/gm  (3*15) 

4 

where  T  *  is  the  temperature  in  units  of  10  degrees . 

Using  p  =  from  (3*3)  and  setting  y 1  =  1.18,  which  is  a  reasonable 

average  (see  Table  II )  we  may  rewrite  (3«13): 


E  =  7.0  X  io11  p"1/^  T'^5 

where  p  is  the  pressure  in  bars.  Since  H  is  proportional  to  E, 
and  ( 3 • 14 )  give 


(3.14) 

(3-12) 


T  ~  ~  p~l*®/(7-l)  =  pj-01 


(3.15) 


Using  7*  -  1  =  0.l8,  which  is  not  far  from  the  average  of  Table  II, 


*The  thermodynamic  relation 

(s£)  =T(l)-p 

T  v 

leads  to  a  relation  between  7-1  and  the  exponents  in  the  relation 
E  =  Ap-X  Ty 


namely. 


x  =  (7  -  1)  (y  -  l) 

Since  we  have  chosen  7  =  1.18  and  y  =  1.5  this  relation  gives  x  =  0.09. 
This  is  in  sufficient  agreement  with  x  =  0.1  as  used  in  (3*13) • 
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(3.16) 


The  numerical  calculations  of  Brode  are  in  good  agreement  with  this  at 
the  relevant  times,  from  about  0.05  to  0.5  second. 

It  will  be  noted  that  (3.16)  was  obtained  without  integration  of 
the  hydrodynamic  equations;  it  follows  simply  from  the  equation  of  state. 
The  weakest  assumptions  are  (l)  the  relation  between  E  and  T,  (3.13) ,  and 
(2)  the  neglect  of  the  constant  in  (3*9).  But  in  any  case,  T  will  be  a 
very  high  power  of  R. 


c.  Mean  Free  Path  and  Radiating  Temperature 


The  emission  of  radiation  from  a  sphere  of  variable  temperature  is 


governed  by  the  absorption  coefficient.  For  visible  light,  the  absorption 

1  ^ 

coefficient  increases  rapidly  with  temperature.  For  any  given  wave 


length,  the  emission  will  then  come  from  a  layer  which  is  one  optical 
mean  free  path  inside  the  hot  material.* 


Actually  the  maximum  emission  comes  from  deeper  inside  the  fireball. 

To  see  this  we  compute  the  emission  normal  to  the  surface,  J  = 

JdrS(R)e  ^dR  where  £.  =  e  emissivity,  and  the  remain - 

cl 

ing  notation  is  as  in  the  text  (below) .  The  integrand  has  a  maximum  at 
R*/l*  =  net  +  Qthv/kT*.  The  optical  depth  at  the  maximum  is  D*  =  D(R*)  = 

□Ct  +  afnv/KT* 

- — — j - which  is  about  2,  rather  than  1  in  the  blue.  By  steepest 

descents  and  with  occasional  use  of  l/nCt  =  0  it  turns  out  that  J  = 

B(v,T*),  greater  by  about  a  factor  2  than  the  J  used  in 

the  text. 
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Assume  that  the  mean  free  path  is 


l 


=  llT 


»-n 


(3.17) 


where  the  exponent  n  ana  the  coefficient  depend  on  the  wave  length. 
Then  the  optical  depth  for  a  given  R  is 


00 


R 


Using  (3-15)  and  (3.17), 


D(R)  = 


1  T 

I 


R 


00 

£[r7  J 


R 


(3.18) 


R 

(nO  -  1)  £(R) 


(3-19) 


A  typical  value  of  n  is  5  [cf.  (5*52)],  and  a  is  10  or  somewhat  larger. 
Thus  to  make 


D(R)  =  1 


(3.20) 


we  need 
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(3.21) 


*(R)  « 

Since  the  interesting  values  of  R  are  of  the  order  of  a  few  hundred 
meters,  the  emitting  layer  will  be  defined  by  the  fact  that  the  mean 
free  path  is  about  5  to  10  meters. 

Equations  (3*18)  to  (3*2l)  hold  for  emission  in  the  exactly  radial 
direction.  At  an  angle  9  with  the  radius  we  get  instead 

*<M)  -  (3.22) 

The  apparent  temperature  of  the  emitting  layer  is  then,  according  to 

(3-17), 

T(R,0)  =  T(R,0)  (cos  S)1/51  (3.23) 

where  T(R,0)  is  the  emission  temperature  for  forward  emission.  The 
intensity  of  emission  is  a  known  (Planck)  function  of  the  wave  length 
and  the  temperature.  Since  the  apparent  temperature  decreases  (though 
slowly)  with  9,  according  to  (3*23)  there  will  be  limb  darkening.  On 
the  many  photographs  of  atomic  explosions,  it  should  be  possible  to 
observe  this  limb  darkening  and  thus  check  the  value  of  n. 

The  relation  (3*17)  needs  to  hold  only  in  the  neighborhood  of  the 
value  (3«2l)  of  l  and  is  therefore  quite  general  as  long  as  &  decreases 
with  increasing  temperature.  The  exponent  n  should  be  determined  at 
constant  pressure.  For  certain  wave  lengths,  especially  in  the  ultra- 
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violet,  (3*17)  is  not  valid;  these  wave  lengths  are  strongly  absorbed  by 
cold  or  cool  air  (Sec.  4c).  For  example,  at  p  =  pQ  and  T  =  2000°,  the 

-j  ^ 

mean  free  path  is  less  than  1  meter  for  all  light  of  hv  >  4.7  ev 
(X  <  0.26  p).  Since  the  emission  of  light  of  such  short  wave  length  from 
such  cool  air  is  negligible,  the  fireball  will  not  emit  such  radiation  at 
all. 

A  detailed  discussion  of  the  absorption  coefficient  in  the  visible 

will  be  given  in  Sec.  4b.  As  can  be  seen  from  the  tables  of  Meyerott 

et  al.1^  and  from  our  Table  VI,  for  a  density  p  =  pQ  the  mean  free  path 

is  of  the  order  of  5  meters  at  about  6000°.  This  corresponds  to  a  pres- 
13 

sure  of  about  25  bars.  For  p  =  O.lp^  the  requisite  mean  free  path  of 
a  few  meters  is  obtained  for  about  10,000°,  with  p  7  bars.  Thus  for  a 
relatively  modest  decrease  in  pressure,  the  effective  temperature  of 
radiation  increases  from  6000  to  10,000°,  corresponding  to  a  very  sub¬ 
stantial  increase  in  radiation  intensity.  This  is  the  mechanism  of  the 
increase  in  radiation  toward  the  second  maximum.  A  more  detailed  dis¬ 
cussion  will  be  given  in  Sec.  5f * 

d.  Energy  Supply 

As  long  as  the  radiating  temperature  is  low,  not  much  energy  will 
be  emitted  as  radiation,  and  this  emission  will  only  slightly  modify  the 
cooling  of  the  material  due  to  adiabatic  expansion.  However,  when  the 
radiating  temperature  increases,  the  radiation  cooling  will  exceed  the 
adiabatic  cooling  to  an  increasing  extent.  It  then  becomes  necessary  to 
supply  energy  from  the  interior  to  the  radiating  surface. 
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Most  of  the  radiation  comes  from  a  thickness  of  one  optical  mean 
free  path  near  the  radius  at  which  ( 3*2l)  is  satisfied.  Let  J  be  the 
radiation  emitted  per  unit  area  per  second  (which  will  be  of  the  order 
of  the  black  body  radiation;  see  below  and  Sec.  4c);  then  the  loss  of 
enthalpy  due  to  radiation,  per  gram  per  second,  will  be 


(3-24) 


Adiabatic  expansion,  according  to  (3«ll),  will  give  an  enthalpy  change 


/  5h\  1  dp 

ydty  p  3t 


(3-23) 


As  long  as  the  shock  is  strong,  i.e.,  as  long  as  p  is  large  compared  to 
the  ambient  pressure  p^,  the  pressure  behaves  as 


p  ~  t 


(3-26) 


where  t  is  the  time  from  the  explosion;  therefore. 


/5h\  _  1. 

\StL,  “  P 


hi  P 
P  t 


(3-27) 


The  radiation  will  be  a  relatively  small  perturbation  as  long  as 
(3-24)  is  smaller  than  (3*27).  This  will  stop  being  the  case  when 


=  1.2 


(3-28) 


Now  using  (3«2l)  this  gives 


(5-29) 


1.2  R 

50  P  t 


=  J 


Because  of  the  steep  dependence  of  temperature  on  R,  ( 3  *  l6 ) ,  the  radi¬ 
ating  surface  R  will  be  close  to  the  shock  front  R  .  How  the  Hugoniot 

s 

relations  state  that  for  y  close  to  1 


R 


(3-30) 


where  p.  is  the  ambient  density  and  p  =  2p,  the  shock  pressure.  Further- 

X  5 

more,  for  the  strong  shock  case. 


=  2.5  $s  (3-31) 

Inserting  into  (3*29), 

0.06  J2  p5//2  =  J  p 11^2  (3-32) 

b  , 

The  black  body  radiation  at  temperature  10  T  is 

Jq  =  5*7  x  1011  T1^  erg/cm2  sec  (3*33) 

Actually,  only  the  radiation  up  to  about  hv^  =  2.75  ev  can  be  emitted  to 
large  distances  because  the  absorption  is  too  great  for  radiation  of 
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higher  frequency  (above,  and  Sec.  i*c).  The  fraction  of  the  black  body 
spectrum  which  can  be  emitted  is  given,  in  sufficient  approximation,  by 


(3.3^) 


where 


0 


(3.35) 


For  T  =  8000°,  we  have  =  4  and  the  effective  emissivity  is  then 

£(4)  =  0.57  (3.36) 

In  (3*34)  we  have  neglected  the  fact  that  the  infrared,  below  about 

l/2  ev,  also  cannot  be  emitted  (Sec.  4a),  and  have  approximated  (eu  -  l) 

in  the  Planck  spectrum  by  e  U;  both  corrections  are  small.  =  8000° 

has  been  chosen  as  a  reasonable  average  temperature  (see  Sec.  5)*  Near 

-1  *5 

this  temperature,  £  varies  about  as  T~  ,  so  that  the  actual  radiation 
to  large  distances  is  about 

J  =  5-7  x  1011  T ^  £(4) 

=  2.3  x  1011  t'2’5  (3.37) 

Solving  (3*32)  for  p,  with  pQ  =  1.29  x  10  J  (normal  air  density)  gives 
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(3-58) 


For  =  pQ,  T*  =0.8,  this  is  14  bars. 

At  higher  temperature  (t'  >  l),  the  ultraviolet  can  be  transported 
as  easily  as  the  visible  although  it  will  not  escape  to  large  distances 
(Sec.  4c).  Then  it  is  reasonable  to  use  the  full  black  body  radiation 
(3*33)  for  the  emission.  Inserting  this  into  (3*32)  gives 


For  p^^  =  pQ,  T*  =  1.0,  this  is  40  bars. 

Thus  for  p  greater  than  14  to  40  bars,  the  radiation  is  only  a 
fraction  of  the  adiabatic  cooling,  for  lower  pressure  radiation  cooling 
is  more  important .  At  the  lower  pressures  then,  energy  must  be  supplied 
from  the  inside  to  maintain  the  radiation.  This  gives  rise  to  a  "cooling 
wave"  moving  inwards  as  will  be  discussed  in  Sec.  5* 

It  is  interesting  that  the  condition  (3*38)  refers  to  the  pressure 
alone.  Neither  the  local  density  nor  the  equation  of  state  enters.  The 
opacity  of  air  enters  only  insofar  as  it  determines  the  radiating  temper¬ 
ature  T*  through  the  condition  (3.21). 

A  more  accurate  expression  for  the  limiting  pressure  will  be  derived 
in  Sec.  5e*  It  will  turn  out  to  be  considerably  lower,  about  5  bars. 
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4.  ABSORPTION  COEFFICIENTS 


a .  Inf raxed 

The  main  absorption  in  the  infrared  is  due  to  free-free  electron 

transitions.  These  are  treated  incorrectly  in  the  paper  by  Meyerott  et 

,  1 6 

al.,  in  which  it  is  assumed  that  such  transitions  occur  only  in  the 
field  of  ions.  At  the  important  temperatures  of  8000°  or  less,  the  degree 

_3 

of  ionization  is  10  or  less.  Therefore,  the  free-free  transitions  in 
the  field  of  neutral  atoms  and  molecules  are  much  more  important  than 
those  in  the  field  of  ions,  even  though  each  individual  atom  contributes 
far  less  than  each  ion. 

The  effectiveness  of  neutrals  in  inducing  free-free  transitions  has 

•j  o 

been  measured  and  interpreted  by  Taylor  and  Fiver  at  the  Avco-Everett 
laboratory.  As  compared  to  one  ion,  the  effectiveness  of  the  most  impor¬ 
tant  neutral  species  is 

N2:  2.2  ±  0.3  X  io"2 

N:  0.9  ±  0.4  X  io"2 

0:  0.2  ±  0.3  x  10~2 

Thus  nitrogen  gives  about  the  same  contribution  whether  molecular  or 
atomic,  and  the  contribution  of  oxygen  is  very  small.  As  a  result,  one 
atom  of  air  is  equivalent  to  about  a  =  0.8  X  io"  ion  (of  unit  charge). 

The  free-free  absorption  coefficient  in  cm"'*'  is  then 

Hff  =  0.87  x  io 5aff,“1/2  2(e")(hv)"5 

=  7.0  t'"i/2  /S-A  2(e")(hv)"3  (4.1) 
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where  e  is  the  number  of  electrons  per  air  atom,  the  quantity  tabulated 
13 

in  Gilmore,  and  hv  is  the  quantum  energy  in  ev.  Table  III  gives  some 
numbers  for  hv  =  1  ev,  four  temperatures  and  three  densities.  At  8000°, 
the  free -free  absorption  is  substantial,  at  lower  temperatures  negligible. 
At  12,000°  the  transitions  are  mostly  in  the  field  of  ions,  i.e., 
Meyerott's  numbers  need  only  slight  correction,  and  the  absorption  is 
large. 

Table  III.  Free -Free  Absorption  Coefficients  (cm~^) 


for  hv  =  1  ev 

T 

p/p0  =  1 

0.1 

0.01 

4,000 

6.8-6 

1.9-7 

4.1-9 

6,000 

5.5-^ 

1.05"5 

r- 

1 

o\ 

• 

1 — 1 

8,000 

4.1-3 

8.7-3 

2.4-6 

12,000 

2.1-1 

1.50"2 

1.12"5 

Note:  for  each  value,  the  power  of  10  is  indicated  by  a  superscript. 

Another  cause  of  absorption  in  the  infrared  is  the  vibrational  bands 

2b  -5  / 

of  NO,  which  have  an  oscillator  strength  of  about  10  and  hv  =  1/4  ev. 
The  resulting  absorption  coefficient  is  about 

=  0.2  (NO)  (4.2) 

where  (NO)  is  the  number  of  NO  molecules  per  air  atom.  This  is  a  few- 
percent  at  p/p0  =  1  and  T  =  4000  to  8000°,  giving  p  =  2  X  10-^  to  10~2 
cm  \  While  this  is  of  the  order  of  magnitude  relevant  for  emission, 
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(3.21),  it  is  small  compared  to  the  free-free  absorption  at  this  low 
frequency,  except  at  T  <  5000°.  Its  main  effect  is  therefore  to  lower 
the  radiating  temperature  in  the  infrared  somewhat. 

Since  the  free -free  transitions  are  the  main  cause  of  infrared 
absorption,  the  electron  density  governs  the  temperature  and  density 
dependence.  For  p/pQ  ^  10 _1  and  T  ^  8000°,  the  main  species  of  positive 
ions  is  NO  .  The  ionization  energy  of  NO  is  9*25  ev;  therefore,  the 
electron  density  is  roughly  proportional  to 


P 


3/2 


exp 


(' 


9*25  ev  \ 
2  kT  J 


(4.3) 


From  6000  to  8000°,  this  gives  a  factor  of  about  10  in  the  electron 
density,  in  accord  with  Gilmore's  tables.  Near  8000°,  we  may  write 
approximately 


Hff  ~  (e“)  ~  p3/2  T7  (4.4) 

1  5 

Writing  p  ~  pT  ,  this  gives 

^ff~P5//2T5  (4.5) 


The  temperature  dependence  of  the  absorption  in  the  visible  will  turn 
out  about  the  same. 


b.  Visible 


In  the  visible  and  for  temperatures  below  about  10,000°,  the  main 


causes  of  absorption  are  the  molecular  bands  of 

Ng  (first  positive)  from  about  1  to  2-l/4  ev 
N*  (first  negative)  "  "  2-l/2  to  3 -3  A  ev 

NO  (P  bands)  "  "  2-l/2  to  6  ev 

and  the  continuous  absorption  due  to  0"  photodetachment  from  about  l-l/2 
ev  up.  The  Ng  second  positive,  and  Og  Schumann -Runge,  although  they  con¬ 
tribute,  are  usually  weaker  than  the  combination  of  N*  first  negative  and 
NO  and  cover  the  same  region  of  the  spectrum  (or  less).  0”  photodetach¬ 
ment  is  mostly  important  at  the  higher  density,  pAq  =  1* 

There  is  commonly  a  'Vindow"  of  low  absorption  below  1  ev,  between 
free-free  absorption  and  Ng(l+)>  and  another  at  about  2~l/2  ev,  between 
Ng  and  the  other  bands  (Table  IV).  The  latter  window  is  filled  in  by  0~ 
photodetachment.  At  higher  temperature,  such  as  12,000°,  the  photoelectric 
effect  on  N  and  0  becomes  important  and  the  molecular  bands  almost  dis¬ 
appear;  the  absorption  coefficient  is  then  almost  uniform  over  the  entire 
spectrum. 

The  table  by  Meyerott  et  al.  should  be  consulted  for  details.  Apart 
from  the  free-free  transitions,  this  table  seems  to  be  in  error  on  the 
Ng(l+)  absorption  at  8000°,  which  should  be  increased  by  a  factor  of  4. 

With  these  corrections,  Table  IV  gives  the  absorption  coefficients  at 
8000°  for  a  few  frequencies  and  densities,  mentioning  in  each  case  only 
the  most  important  species  of  absorbers.  For  pAq  =  0.1  and  0.01,  the 
contribution  of  NO  and  0g  in  the  wave  length  region  considered  is  negli¬ 
gible.  The  windows  at  1  and  above  2  ev  are  noticeable  in  the  table;  beyond 
3  ev  there  is  a  rapid  increase  in  absorption. 
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Table  IV.  Contributions  to  Absorption  Coefficients  at  8000° 


hv  (ev) 

5/8 

1-1/8 

1-5/8 

2-1/8 

2-5/8 

3-1/8 

p/p0  -  i 

ff 

1.7"2 

0.3"2 

0.10-2 

0.04"2 

0.02"2 

0.01"2 

l.l-2 

2.1"2 

0.84"2 

0.04"2 

0.56"2 

O' 

0.21-2 

0.28“2 

0.31-2 

0.32"2 

NO,  Ng,  02 

_p 

0.31 

2.2'2 

Total 

1.7"2 

1.4'2 

2.4"2 

1.16'2 

0.68"2 

3.1-2 

p/p0  =  0*1 


ff 

0.36"5 

0.06 “5 

0.02-5 

O.Ol"5 

N2 

0.6l"5 

1.12-5 

1 

—4 

• 

0 

0.02-5 

0.30-5 

0“ 

o.05“5 

0.0 6-5 

0.07~5 

0.07"5 

A 

0.55-5 

4.2"5 

Total 

0.36"5 

0.67-5 

1.2-5 

0.53"5 

0.64"5 

4.6'5 

p/p0  =  0.01 


ff 

l.o"5 

0.17~5 

0.05-5 

0.02  "5 

0.0l"5 

O.Ol"5 

N2 

1.39-5 

2.6-5 

1.03-5 

0.05-5 

0" 

0.13-5 

0.17-5 

0.2  "5 

0.2"5 

K2 

0.03-5 

4.4-5 

34-5 

Total 

l.O"5 

1.6"5 

L C\ 

I 

CO 

OJ 

1.25-5 

4.7-5 

SO'5 
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The  concentration  of  all  the  absorbing  species  depends  strongly  on 
temperature:  The  Ng  first  positive  absorption  starts  from  the  electronic 
level  A,  which  has  an  excitation  energy  of  5*7  ev;  NO  itself  requires 
high  temperature  for  its  formation;  and  Ng  and  0~  are  ions  whose  concen¬ 
tration  behaves  like  the  electron  concentration,  discussed  above  in 
Sec.  4a.  Therefore,  the  absorption  depends  strongly  on  temperature .  Ng 
first  positive  is  the  most  important  absorber,  and  its  temperature  depend¬ 
ence  (relative  to  the  ground  state  of  Ng)  is  about 


The  concentration  of  Ng  -up  to  6000°  is  nearly  independent  of  density, 
but  at  8000°  it  is  about  proportional  to  p^2  (per  air  atom)  so  that  near 
8000° 


~  p5//2  T8  (It. 7) 

2 

which  is  nearly  the  same  dependence  as  derived  in  (4.4)  for  the  free-free 
transitions.  In  terms  of  pressure  we  get 

~  P3^2  T6  (4.8) 

iN2 

Table  V  gives  the  absorption  coefficients  for  four  temperatures, 
three  densities,  and  six  wave  lengths.  In  the  visible  (l-l/8  to  2-5/8  ev) 
the  strong  increase  in  absorption  with  temperature  is  evident,  a  factor 
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Table  V.  Absorption  Coefficients  in  Visible 


CO 

O 

KP 

1 

CM 

1 

tP- 

1 

O 

UP  -4 

CM 

1 

OJ 

CO 

1 

VO 

UP 

KP 

l 

H 

H 

O 

H 

O 

1  l 

O 

1 

1 

CP 

o\ 

pH 

• 

pH 

VO  vo 
cm  -4 

• 

O 

UP 

-4 

• 

KP 

CO 

1— 1 

KP 

-4 

rH 

CP 

CM 

KP 

O 

CO 

VO 

KP 

1 

VO 

OJ 

1 

CO 

tp- 

UP 

1 

-4 

OJ 

up 

CO 

VO 

UP 

1 

KP 

1 

00 

UP 

vo 

vo 

1 

CO 

-4 

co 

1 

CO 

t- 

vo 

1 

OJ 

• 

• 

UP 

• 

• 

• 

vo 

• 

* 

• 

CM 

UP 

H 

0 

KP 

-4 

vo 

0 

H 

CM 

*4 

0 

KP 

1 

0J 

t- 

CM 

CO 

UP 

KP 

CO 

vo 

1 

UP 

-4 

1 

1 

vo 

1 

1 

14- 

H 

l 

UP 

KP 

UP 

© 

UP 

l 

KP 

5$ 

0 

OP 

1 

CO 

UP 

0J 

1 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

CM 

H 

H 

H 

pH 

IP- 

UP 

O 

0 

-4 

pH 

0 

CO 

VO 

KP 

CM 

1 

UP 

*KP 

H 

-4 

1 

• 

• 

• 

H 

UP 

CM 

CM 

KP 

i 


ip- 

UP 

CM 

CO 

vo 

UP 

1 

1 

1 

1 

I 

1 

1 

O 

CM 

14- 

-4 

UP 

O 

KP 

CO 

KP 

• 

• 

l 

• 

• 

• 

• 

• 

-4 

KP 

OJ 

pH 

-4 

O 

CM 

H 1 

pH 

pH 

pH 

CO 

1 

O 

OJ 

1 

tP“ 

1 

UP 

l 

-4 

1 

CM 

CO 

1 

VO 

1 

UP 

1 

CO 

pH 

tr\ 

pH 

-4 

CM 

vo 

t4- 

OJ 

KP 

IP- 

vo 

-4 

1 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

rH 

\ — 1 

pH 

-4 

tp- 

vo. 

pH 

KP 

UP 

pH 

O 

CO 

IPs 

II 

li 


KP 

I 

KP 

CM 


t— 

fl 

UP 

1 

-4 

1 

CM 

1 

CO 

1 

UP 

1 

O 

KP 

1 

CO 

KP 

vo 

pH 

IP- 

00 

VO 

• 

• 

• 

• 

• 

m 

• 

• 

KP 

VO 

1— 1 

0 

rH 

-4 

8 

O 

O 

8 

O 

O 

8 

O 

O 

8 

O 

O 

O 

Q 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

-4* 

vcT 

co' 

vcT 

co" 

£ 

.4* 

vo' 

co" 

pH 

O 


O 

o 
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of  about  200  from  4000  to  6000°,  about  10  from  6000  to  8000°  except  for 

the  lowest  density,  and  10  to  50  from  8000  to  12,000°.  In  most  cases, 

the  dependence  on  wave  length  is  slight  until  the  rapid  increase  of 

absorption  in  the  ultraviolet  which  sets  in  at  about  2.5  ev  for  4000°, 

5  ev  at  6000  and  8000°,  and  not  at  all  at  12,000°.  The  especially  strong 

increase  for  4000°  is  due  to  the  Schumann -Runge  bands,  which  are  not  very 

*,t  ,  ,  . 

sensitive  to  temperature  (see  Sec.  4c). 

In  Sec.  5d  we  shall  need  the  mean  free  path  in  gm/cm  ,  suitably 
averaged  over  the  "transparent"  region.  From  Table  V  it  appears  that  a 
reasonable  estimate  of  this  region  is  from  hv  =  l/2  to  2-5/4  ev.  We 
have  averaged  p/p  as  calculated  from  Table  V,  with  the  weighting  factor 


3  -u 
u  e  , 


k  X  8000° 
hv 


(4.9) 


which  does  not  vary  much  (from  0.84  to  1.33)  between  l-l/8  and  2-5/8  ev. 
The  result  is  given  in  Table  VI. 


*After  completion  of  this  paper,  I  received  new  absorption  coefficients 
by  F.  R.  Gilmore,  cf.  ref.  23*  Aside  from  including  the  absorption  by 
NO2  (see  Sec.  4a),  Gilmore  includes  the  free-free  absorption  in  the 
field  of  neutrals,  based  on  ref.  l8,  and  also  new  data  on  f  numbers  for 
the  important  bands. *4  The  most  important  change  is  a  reduction  in  the 
f  number  of  the  W2  first  positive  system  from  0.02  to  0.0028,  which  will 
substantially  reduce  the  absorption  in  the  visible.  Unfortunately,  this 
will  further  raise  the  theoretical  radiating  temperature  (Sec.  5d) , 
which  is  already  higher  than  observed. 

^Some  recent  Avco-Everett  experiments  may  indicate  that  the  free-free 
electron  transitions  in  the  field  of  the  N  atom  are  enhanced  in  the 
visible  as  compared  to  the  infrared.  This  may  compensate  to  some  extent 
for  the  reduction  of  the  Ng  first  positive  bands  (see  footnote  above). 
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Table  VI.  Average  Mean  Free  Paths  (gm/cm  ) 


T 

H 

11 

O 

a 

0.1 

0.01 

6,000 

0.86 

1.2fc 

3.15 

8,000 

0.105 

0.20 

0.76 

12,000 

_2 

1.25 

1.82': 

A  fair  approximation  to  Table  VI  is 


c.  Ultraviolet 

The  absorption  in  the  ultraviolet  is  generally  high  at  all  tempera¬ 
tures.  At  low  temperatures,  the  absorption  is  mainly  due  to  the  Schumann  - 
Runge  bands;  at  higher  temperature  (8000°)  these  are  replaced  by  NO  3  and 
y,  and  at  high  temperature  (12,000°)  by  photoelectric  absorption  in  0~, 


h-0 


N,  and  0.  The  Schumann -Runge  bands  start  from  the  electronic  ground  state 
of  02,  hence  are  available  at  low  T;  at  higher  T,  the  spectral  region  of 
strong  absorption  spreads  due  to  excitation  of  vibrational  states;  but  at 
still  higher  T,  oxygen  dissociates  and  therefore  the  bands  die  out  (at 
8000°,  they  contribute  less  than  10$  of  the  absorption).  The  photo¬ 
electric  absorption  in  N  and  0  requires  not  only  the  presence  of  these 

‘X* 

atoms  but  also  their  electronic  excitation,  and  therefore  does  not  become 
important  until  about  10,000°. 

Table  VII  gives,  for  density  p  =  O.lp^,  the  spectral  regions  of 

_3 

strong  absorption.  In  accord  with  Sec.  3b,  we  define  this  by  p  >  10 

-2 

(mean  free  path  less  than  10  meters)  or  |i  >  10**  (£  <  1  meter).  The 

table  shows  that  strong  absorption  covers  a  particularly  wide  spectral 
region  at  4000°,  shrinking  substantially  at  6000°.  Very  strong  absorption, 
\i  >  10 cm-^,  occurs  in  quite  a  large  spectral  region  for  T  =  ^000°  but 
shrinks  to  practically  nothing  at  6000°  and  to  nothing  at  8000°.  At 
12,000°,  very  strong  absorption  occurs  again  but  is  now  in  the  visible. 

The  strong  absorption  in  the  ultraviolet  (hv  >  3«5  ev)  means  first 
of  all  that  the  UV  is  not  emitted  to  large  distances  and  can  therefore 
not  be  observed;  e.g.,  at  T  =  4000°  and  hv  =  3*5  ev,  we  have 


hV  m 

u  ”  a  ‘  10 


*We  should  also  add  the  photoelectric  absorption  from  excited  states  of 
NO,  which  should  be  especially  noticeable  at  8000°. 


The  fraction  of  the  Planck  spectrum  beyond  u  =  10  is  only  about  1 $>t  so 
that  emission  of  these  frequencies  is  negligible. 

Table  VII.  Ultraviolet  Absorption:  Spectral  Regions  (hv  in  ev)  with 

Large  Absorption  as  a  Function  of  Temperature  for  p  =  0.1pQ 

T  \i  >  10  "5  cm"1  P  >  10  "2  cm-1 


2,000 

4.7  -  7.2 

5.5  -  7.2 

3,000 

3.9  -  7.2 

4.7  -  7.2 

4,000 

3*5  -  7.2 

4.6  -  7*2 

6,000 

4.0  -  7.1 

5.8  -  6.0 

8,000 

2.7  -  6.3 

None 

12,000 

All 

2.7  -  3-5 

The  ultraviolet  can,  however,  be  transported  quite  easily  at  8000° 
and  even  more  easily  at  12,000°  if  there  is  a  temperature  gradient.  Such 
a  gradient  is  always  available,  whether  we  have  adiabatic  conditions 
(Secs.  3b,  5f)  or  a  strong  cooling  wave  (Sec.  5d) .  Therefore,  there  will 
be  a  flow  of  ultraviolet  radiation  at  the  radiating  temperature,  defined 
in  Sec.  5,  which  will  be  shown  (Secs.  5d,  5f )  to  be  about  10,000°  or 
slightly  less.  To  calculate  this  flow,  we  should  determine  the  tempera¬ 
ture  gradient  from  considerations  such  as  Sec.  5d  or  5f,  and  then  insert 
this  into  the  radiation  flow  equation.  This  is  similar  to  (5.3)  except 
that  only  the  ultraviolet  contribution  to  the  flow  should  be  taken  into 


account 


When  this  is  done  in  a  case  of  constant  ( frequency -independent ) 
absorption  coefficient,  the  ultraviolet  transport  will  he  related  to  the 
visible  radiation  transport  as  the  respective  intensities  in  the  Planck 
spectrum.  This  condition  seems  to  be  nearly  fulfilled  at  12,000°.  At 
8000°,  the  absorption  in  the  near  ultraviolet  (2.75  to  4.2  ev)  is  about 
three  times  that  in  the  visible;  then  the  UV  transport  will  be  one -third 
of  that  corresponding  to  the  Planck  intensity.  Since  the  radiating 
temperature  near  the  second  radiation  maximum  is  between  8000  to  10,000°, 
the  actual  UV  transport  will  be  between  one -third  and  the  full  Planck 
value,  relative  to  the  visible  radiation.  According  to  (3*56),  the  UV 
contains  about  43$  of  the  Planck  intensity  at  8000°;  hence  the  total 
radiation  transport  at  this  temperature  is  about 

57  +  1/3  x  45  =  71$ 

of  the  black  body  radiation.  At  12,000°  we  get  the  full  black  body  value. 
For  simplicity  we  have  assumed  the  full  black  body  radiation  in  Sec.  5, 
even  though  the  UV  is  not  emitted  to  large  distances.  But  this  problem 
could,  and  should,  be  treated  more  accurately. 

Having  discussed  the  influence  of  the  UV  on  the  total  radiation  flow, 
we  now  examine  what  happens  to  the  UV  radiation  after  it  has  gone  through 
the  "radiating  layer,"  i.e.,  the  layer  which  emits  the  visible  light  to 
large  distances.  The  very  near  ultraviolet,  2.75  to  3*5  ev,  will  be 
partially  absorbed  at  4000  to  6000°,  especially  if  the  layer  of  matter 
at  these  intermediate  temperatures  becomes  thick,  0.3  to  0.5  gm/cm  or  so. 
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The  UV  beyond  3.5  ev  will  be  strongly  absorbed  at  4000°.  Thus  the  layers 
of  air  at  intermediate  temperatures  get  additional  heat  which  counteracts 
and  may  even  exceed  the  adiabatic  cooling.  This  will  tend  to  increase 
the  thickness  of  the  medium -temperature  layer.  This  in  turn  will  (mod¬ 
erately)  lower  the  radiating  temperature,  but  three-dimensional  effects 
act  the  opposite  way  (Sec.  6). 

5.  THE  COOLING  WAVE 
a.  Theory  of  Zel'dovich  et  al. 

Zel'dovich,  Kompaneets,  and  Raizer  ^  (quoted  as  Z)  have  considered 
the  loss  of  radiation  by  hot  material  when  the  absorption  coefficient 
for  the  radiation  increases  monotonically  with  temperature.  They  have 
shown  that  in  this  case  a  cooling  wave  proceeds  into  the  hot  material 
from  the  outside.  This  is  to  say,  the  cool  temperature  outside  gradually 
eats  its  way  into  the  hot  material,  while  the  material  in  the  center 
remains  unaffected  and  merely  expands  adiabatically. 

For  simplicity,  Zel'dovich  et  al.  consider  a  one -dimensional  case. 
They  further  assume  that  the  specific  heat  is  constant  and  express  their 
theory  in  terms  of  the  temperature.  This  is  not  necessary;  we  shall 
merely  assume  that  both  the  enthalpy  H  and  the  absorption  coefficient 
for  radiation  are  arbitrary  but  monotonically  increasing  functions  of  the 
temperature.  Like  Z,  we  shall  assume,  in  this  subsection  only,  that  the 
radiation  transport  can  be  described  by  an  opacity  (Rosseland  mean)  rather 
than  considering  each  wave  length  separately. 
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The  fundamental  statement  of  Z  is  that  the  cooling  wave  keeps  its 
shape,  i.e. ,  that  the  enthalpy  (and  other  functions  of  the  temperature) 
is  given  by 

H  =  H(x  +  ut)  (5.1) 


Here  t  is  the  time,  x  is  the  Lagrangian  coordinate,  and  u  is  the 
Lagrangian  velocity  of  the  cooling  wave.  We  have  written  x  +  ut  so  that 
the  cooling  wave  proceeds  towards  smaller  x,  i.e.,  inwards.  H  is,  of 
course,  a  decreasing  function  of  x  +  ut.  The  Lagrangian  coordinate  is 
best  measured  in  gm/cm  and  is  defined  by 


(5-2) 


where  X  is  the  geometrical  (Eulerian)  coordinate.  For  given  pressure  p, 
the  density  p  is  a  function  of  H  so  that  X(x)  can  be  calculated  from 
(5*2).  The  Lagrangian  velocity  u,  measured  in  gm/cm  sec,  is  a  constant. 

For  any  given  H,  we  know  the  temperature  T,  hence  the  opacity  K  and 
the  radiation  flow 


4  a  d(T4 5) 

5  K(T)  ~<5x 


(5-5) 


where  a  is  the  Stefan-Boltzmann  constant 


c  2 

=  5*7  x  10~^  erg/cm  sec  deg 
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The  equation  of  continuity  is 


SH  _  dJ 
dt  §x 


(5.4) 


The  energy  in  radiation  has  been  neglected,  which  is  justified  in 
all  practical  cases.  Using  (5»l),  (5*4)  can  be  integrated  over  x  to 
give 


uH  +  J  =  C  (5.5) 

where  C  is  a  constant.  This  is  the  fundamental  result  of  Z. 

If  the  opacity  increases  monotonically  with  H,  then  in  the  interior 
J  will  be  very  nearly  zero,  and  therefore 

C  =  uHq  (5.6) 


where  HQ  is  the  enthalpy  In  the  undisturbed  interior,  hot  region. 
Equation  (5*5)  becomes 


J  =  u(Hq  -  H) 
and  using  (5*5) 


(5.7) 


dx  = 


4  d(T4) 

5  a  K(T,p)(H0  -  H) 


(5.8) 


which  can  be  integrated  to  give  x(T),  since  H(T)  and  K(T,p)  are  known 
functions.  We  have  put  in  evidence  the  fact  that  K  depends  on  pressure 
in  addition  to  T.  Over  most  of  the  range  of  T,  K(T)  is  the  most  rapidly 
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varying  (increasing)  function  of  T  and  the  variation  of  -  H  is  less 

important;  therefore,  T(x)  becomes  steeper  as  T  increases;  but  when  H 

gets  very  close  to  H^,  the  most  rapidly  varying  function  in  (5*8)  is 

OCX 

Hq  -  H,  and  H  approaches  HQ  exponentially  as  e  for  small  x.  The  qual¬ 
itative  behavior  of  T(x)  is  shown  in  Fig.  1,  in  accord  with  Z.  To  obtain 

X 

this  shape  it  is  essential  that  K(T)  increase  much  faster  than  T  . 


H0 


Fig.  1.  Temperature  distribution  in  cooling  wave. 

On  the  outside,  we  finally  come  to  a  point  x^  where  only  one  optical 
mean  free  path  is  outside  x^.  From  this  point  we  get  black  body  emission, 
i.e.,  (5*3)  is  replaced  by 

J(Xi)  =  Ji  =  aT^  (5-9) 


Using  this  in  (5«7)  we  find 


u 


(5.10) 


v^r 

To  determine  u  we  therefore  have  to  proceed  as  follows: 

1.  Find  the  temperature  T^  at  which  the  opacity  K(T^)  is  such  that  there 
is  one  optical  mean  free  path  outside  i.e., 

T1 

J  K(T)  dx(T)  =  1  (5.H) 

0 

For  this  purpose  we  must  know,  of  course,  the  temperature  distribution 
T(x)  for  T  <  T^. 

2.  Determine  J(T^)  from  (5«9)  and  H(T^)  from  the  equation  of  state. 

Knowing  the  internal  enthalpy  HQ  then  gives  u  from  (5.10).  Note 
that  u  is  the  Lagrangian  velocity  of  the  cooling  wave.  It  has  the  cor¬ 
rect  dimension. 

To  solve  problem  1,  Z  assume  that  the  material  which  has  gone  through 
the  cooling  wave  will  expand  adiabatically.  We  shall  find  that  this  is  a 
reasonable  assumption  in  most  conditions  (Sec.  5d)  but  that  at  early  times 
(Sec.  5f)  and  in  certain  late  stages  other  considerations  apply  (Sec.  6b). 

b.  Inside  Structure  of  Fireball,  Blocking  Layer 

In  early  stages  (Stage  B  I),  just  after  the  shock  wave  is  formed, 
the  isothermal  sphere  expands ,  by  radiation  diffusion,  into  the  material 
which  has  been  heated  by  shock.  This  process,  which  will  be  treated  in 
a  subsequent  report,  depends  on  the  temperature  and  temperature  gradient 
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in  the  isothermal  sphere.  For  its  occurrence  it  is  Important  that  the 
opacity  actually  decreases  with  increasing  T  at  higher  temperature.  Many 

calculations  of  K(T)  in  this  temperature  region  have  been  made.  Curves 

21 

have  been  compiled,  especially  by  Gilmore,  and  revised  as  more  infor¬ 
mation  has  become  available.  The  most  recent  and  most  extensive  calcu- 
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lation,  to  my  knowledge,  is  by  Stuart  and  Pyatt. 

All  calculations  agree  that  T^/k(T) ,  which  is  the  important  quantity 
according  to  (5*3)>  has  a  pronounced  minimum  at  about  T*  =  2,  T  =  20,000°. 
(The  temperature  of  the  minimum  increases  slightly  with  increasing  den¬ 
sity.)  Molecules  are  no  longer  present  at  these  temperatures,  and 
absorption  is  mainly  by  bound-free  electron  transitions  in  atoms  and 
atomic  ions,  with  some  contribution  from  broadened  atomic  lines  ( bound - 
bound  transitions)  whose  calculation  is  the  most  difficult.  Radiation 
transport,  then,  is  most  difficult  around  T^  =  20,000°,  and  temperatures 
around  T^  constitute  an  effective  blocking  layer  for  radiation. 

Until  about  0.5  second,  the  central  temperature  of  a  1-megaton  sea 

level  explosion  is  greater  than  T^,  according  to  the  calculations  of 

8 

Brode.  Radiation  flow  can  then  be  considered  as  taking  place  separately 
in  an  interior  and  an  exterior  region.  The  interior  flow  is  determined 
by  the  central  temperature  T  ,  and  this  flow  generally  decreases  with 
time  because  T  decreases  and  with  it  the  quantity  T^/k.  The  exterior 
flow  is  dominated  by  the  cooling  wave  and  increases  with  time  because 
the  decrease  of  density  causes  a  decrease  of  opacity  for  any  given  T. 

The  two  flow  regions  are  separated  by  a  blocking  layer  in  which  the 


temperature  is  around  T^,  and  in  which  the  temperature  distribution  is 
essentially  that  originally  established  by  shock  and  subsequent  adiabatic 
expansion.  It  is  to  this  condition  that  the  temperature  distribution 
(5*l6),  T  ~  R-"^,  refers. 

The  radiation  flow  through  the  blocking  layer  is 


J  =  .  i  aT4  ^  4  d  lQS  T 
d  3  b  Rb  4  d  log  R 


(5.12) 


The  mean  free  path  for  radiation  in  the  blocking  layer  (at  18,000°) 
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according  to  Gilmore  is  about 


cm 


(5.13) 


The  radius  R^  at  which  this  temperature  occurs,  for  a  1-megaton  explosion 
at  sea  level,  is  R^  =  300  to  1+00  meters.  Using  T  ~  R~^,  this  gives 


7  x  109 


(5.1M 


This  is  equivalent  to  black  body  emission  at  the  effective  temperature 


Teff,b  *  5500 


(y 


-0.325 


(5.15) 


Using  the  equation  of  state  at  18,000°, 
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(5.16) 


where  p  is  in  bars.  Equation  (5*15)  becomes 

Teff,b  =  20 >000  p’0,36  (5.17) 

The  "blocking  layer"  will  deserve  its  name  only  if  the  radiating  temper¬ 
ature  is  higher  than  T  ,  . 

0X1  y  ID 

c.  Velocity  of  Cooling  Wave 

The  velocity  of  the  cooling  wave  is  given  by  (5.10),  where  is 
the  enthalpy  at  the  point  to  which  the  cooling  wave  has  proceeded.  If 
inside  this  point  there  is  a  noticeable  flow  of  radiation,  J ,  (5*10) 
should  be  generalized  to 


u 


(5.18) 


Usually,  the  main  dependence  on  the  internal  conditions  is  through  HQ, 
the  effect  of  JQ  being  less  important. 

As  Hq  increases,  i.e.,  as  the  wave  progresses  more  into  the  interior, 
the  velocity  of  the  wave  will  decrease.  The  limit  will  be  reached  when 
the  cooling  wave  penetrates  the  isothermal  sphere;  then  is  the  enthalpy 
in  that  sphere  and  JQ  =  0.  We  now  use  the  black  body  formula  (5*55)  for 
the  emission  of  radiation  at  the  radiating  surface  T^,  (3*1*0  for  the 
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internal  energy  in  the  isothermal  sphere  of  temperature  Tq  and  y'  -  1.15 
in  that  sphere;  then  (5-l8)  becomes 


u 


5.7  x  1011  T.^ 


°  8.2  x  1011  p“1,/5  -  T*^5) 

c  1 


=  0.70  p 


1/9 


rp'5/3  _  rp'5/3 

C  1 


sec 


(5.19) 


Typically,  p  =»  5,  =  1,  «  3;  then  u  =  0.15  gm/cm2  sec.  The  expres¬ 

sion  (3.33)  includes  all  the  black  body  radiation.  If  only  the  radiation 
actually  emitted  to  large  distances  is  to  be  included  (which  is  reasonable 
at  lower  temperatures,  T^  <  0.8)  (3*37)  should  be  used;  then  u  will  be 
smaller,  0.1  gm/cm  sec  or  less. 

Before  the  cooling  wave  reaches  the  isothermal  sphere,  HQ  is  smaller. 
This  is  partly  compensated  by  the  fact  that  JQ  >  0.  An  interesting  inter¬ 
mediate  state  is  when  the  cooling  wave  has  just  reached  the  blocking 
layer.  Then,  using  T^  =  1.8  and  (5.1k), 


(5.20) 


=  0.70 


16 


-l.hb 

P 


7^73 

1 


(5.21) 


Clearly  this  makes  sense  only  if  the  subtracted  term  in  the  numerator  is 
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smaller  than  the  first  term,  i.e.,  when  the  pressure  is  sufficiently 
high.  As  the  pressure  decreases  belcw  about  10  bars,  the  blocking  layer 


"opens  up"  and  ceases  to  block  the  flow  of  radiation. 

In  the  very  beginning,  when  the  cooling  wave  just  starts,  the  top 
temperature  of  the  cooling  wave,  TQ,  is  close  to  the  radiating  tempera¬ 
ture  T  .  Then  (5*18)  becomes 

"i  -  -  (§1  (5'22) 

1 

If  we  use  (5.12)  for  J,  assume  d  log  T/d  log  R  and  R  to  be  constant,  and 
use  (3.14),  then 


For  further  discussion,  see  Sec.  5f* 

As  T  increases,  u  decreases  from  (5*23)  via  (5*2l)  to  (5*19)*  After 
the  cooling  wave  has  penetrated  to  the  isothermal  sphere,  u  is  apt  to 
increase  again  because  Tc  in  the  denominator  of  (5*19)  will  decrease  due 
to  adiabatic  expansion  of  the  isothermal  sphere.  Thus  the  velocity  u  is 
apt  to  be  a  minimum  when  the  cooling  wave  has  just  reached  the  isothermal 
sphere . 

The  variation  of  u  with  time  is  not  very  great.  Likewise,  the  shape 
of  the  cooling  wave  changes  only  slowly  with  time.  The  shape  is  obtained 
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by  integrating  (5*8);  it  depends  on  time  because  K  is  a  (rather  slowly 
variable)  function  of  the  pressure,  and  Hq  is  a  (slow)  function  of  time. 
This  justifies  approximately  the  basic  assumption  (5.1)  of  Z:  While  the 
cooling  wave  does  not  preserve  its  shape  exactly,  it  does  so  approximately. 

The  picture  is  then  that  at  any  given  time  t  there  are  up  to  five 
regions  behind  the  shock.  In  Fig.  2  the  temperature  is  plotted  schemat¬ 
ically  against  the  Lagrange  coordinate  r.  Starting  from  the  center, 
there  is  first  the  isothermal  sphere  (i  in  Fig.  2). 


Fig.  2.  Schematic  temperature  distribution.  I,  isothermal 
sphere;  III,  cooling  wave;  VI,  undisturbed  air;  E, 
shock  front.  II,  IV,  and  V  are  expanding  adiabat- 
ically. 


This  may  be  followed  by  a  region  II  in  which  the  temperature  distribu¬ 
tion  is  essentially  that  established  by  adiabatic  expansion  behind  the 
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shock,  Eq.  (3. 15).  Next  comes  the  cooling  wave  III  in  which  the  temper¬ 
ature  falls  more  steeply,  according  to  (5*8).  (At  late  times,  region  II 
is  wiped  out  and  III  follows  immediately  upon  I.)  Region  IV  includes 
the  material  which  has  gone  through  the  cooling  wave,  and  now  cools  adi- 
abatically;  hence  the  temperature  falls  slowly  with  r  (Sec.  5d).  D  is 
the  material  point  from  which  the  cooling  wave  started  originally;  region 
V,  outside  that  point,  is  also  expanding  adiabatically,  but  from  shock 
conditions;  thus  it  is  the  continuation  of  region  II.  Finally  region  VI 
is  the  air  not  yet  shocked.  As  time  goes  on,  the  cooling  wave  III  moves 
inward,  wiping  out  region  II  and  then  eating  into  region  I.  Region  IV 
accordingly  grows  toward  the  inside,  but  its  outer  end  D  stays  fixed. 
Region  V  expands  into  VI  by  shock. 

We  note  once  more  that  u  is  the  velocity  in  Lagrange  coordinates, 
and  in  gm/cm2  sec.  The  problem  is  made  somewhat  more  complicated  by  the 
three  dimensions  and  the  adiabatic  expansion,  cf.  Sec.  6,  but  the  prin¬ 
cipal  features  remain  the  same. 


Adiabatic  Expansion  after  Cooling.  Radiating 


Derat ure 


When  a  given  material  element  has  gone  through  the  cooling  wave,  it 


is  left  at  the  radiating  temperature  T^.  Thereafter,  it  will  expand 
adiabatically.  Equation  (3.25)  shows  that  for  adiabatic  expansion 


(  9h\  1  dp 

Wadi  =  ?  * 


(5.24) 


or,  using  (3*10) 
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9  log  II  _  7'  -  1  9  log  p 
<5t  71  9t 


(5.25) 


As  long  as  the  shock  is  strong,  (3*26)  holds  and  therefore 


5  log  H  _  _  .  p  7  -  1 
9  log  t  ~  y ' 


(5.26 


Gilmore  * s  J  Tables  11  and  15  show  that  for  T  =  5000  to  8000°,  and  p  = 
1  to  10  bars,  y*  varies  from  about  1.13  to  1.20,  so 


9  log  H 
9  log  t 


=  -  0.1'!-  to  -  0.20 


(5.27) 


At  late  times,  p  no  longer  decreases  as  fast  as  t"  *  ,  so  H  also  decreases 
more  slowly  but  (5*25)  remains  valid. 

We  use  now  the  approximate  equation  of  state  (3*li0  (together  with 
H  =  7' E)  and  find 


v  ith 


T'~Pl/l5  h3/5 


_L  ,  3  y’-i 

15  5  ~y~ 


(5.28) 


2  7  -  0.9 

5  p 


0.136  to  0.167 


(5.29) 
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using  (5*25)  and  7*  =  1.15  to  1.20.  A  given  material  element  which  went 

through  the  cooling  wave  at  temperature  T '  and  pressure  p  will  now  (at 

m  *  rm 

pressure  p)  have  the  temperature 


(5-30) 


If  we  assume  T  to  be  independent  of  time,  then  at  a  given  time  in  the 
adiabatic  region  (IV  of  Fig.  2) 


3  log  I'  _  paiogt»_Sa  lc« 
dx  5x  u  5t 

m 


(5*31) 


1.23  _  0^8  7 1  -0.9 

ut  "  ut  71 
m  m  ' 


(5-32) 


where  tffi  is  the  time  at  which  the  material  element  x  was  radiating.  The 
last  step  in  (5«3l)  assumes  that  the  velocity  of  the  radiation  wave  is 

constant  but  is  valid  whatever  the  time  dependence  of  the  pressure;  (5.32) 

-1*2 

assumes  p~t  ,  i.e.,  strong  shock  conditions.  It  should  be  noted 
that  t  in  the  relation  p  ~  t-1*2,  and  therefore  t  in  (5.32),  is  the  total 
time  since  the  nuclear  explosion. 

The  velocity  u  of  the  cooling  wave  was  calculated  in  Sec.  pc,  (5«l8), 
(5.19),  etc.  The  condition  for  the  radiating  surface,  now  as  in  Sec.  3c, 
is  that  there  be  one  optical  mean  free  path  outside  it  for  "visible"  light 
as  defined  in  Sec.  4b.  Since  we  now  use  material  coordinates  in  gm/cm2, 
we  should  use  mean  free  paths  in  the  same  unit.  Table  VI  and  (4.10)  and 
(4.12)  give  the  required  information.  We  write  (4.12)  in  the  form 
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£  -  Ap-1/3 

n 


(5.33) 


2 

with  A  =  0.10  gm/cm  and  n  =  6.5. 

The  optical  depth  is,  according  to  (5*32)  and  (5 <>33) 


Setting  D  =  1,  7*  =  1.15,  A  =  0.10,  and  n  =  6.5  gives 


rn’6‘5  _  0.11  gm/cm2  -l/5 
1  ut  P 


(5.36) 


This  is  an  explicit  expression  for  'the  radiating  temperature  in  terms  of 
the  velocity  of  the  cooling  wave.  The  latter  depends  in  turn  on  the 
radiating  temperature,  increasing  with  T^  ,  so  that  T  occurs  altogether 
in  the  10.5  power  and  thus  can  he  determined  very  accurately. 

Equation  (5«3 6)  can  he  further  reduced  hy  using  the  relation  between 
pressure  and  time  which  is,  for  a  strong  shock,  approximately 


(5.37) 
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where  t  is  in  seconds,  p  in  bars,  and  Y  in  megatons,  and  is  the  den¬ 
sity  of  the  ambient,  undisturbed  air.  Then  (5.36)  becomes 


(5.38) 


The  right  hand  side  of  (5*38)  gives  the  complete  dependence  on  Y  and 
since,  in  deriving  (5*38),  we  have  only  used  the  opacity  law  (5*33)  and 
the  adiabatic  cooling  of  air,  (5*32),  both  of  which  are  independent  of 
the  explosive  yield  Y  and  of  the  ambient  density  of  the  air. 

Now  insert  u  from  (5*19);  then  we  obtain 


m'10.5 

1 _ 

t*5/3  _  t'5/3 
c  1 


=  0.16  Y 


-1/3 


p7/l8 


(5-39) 


This  equation  gives  the  radiating  temperature  in  terms  of  the  central 
temperature  T^  and  of  the  quantities  on  the  right  hand  side.  The  radi¬ 
ating  temperature  is  proportional  to  a  low  power  of  the  central  temper¬ 
ature  (about  the  l/6  power);  thus  as  the  inside  cools,  the  radiation 
decreases  (see  Sec.  5e  for  details).  It  also  decreases  slowly  with  time 
due  to  the  pressure  factor  on  the  right  hand  side,  T^  ~  pl/27 ^  por  gjj_ven 
p  and  Tc,  the  radiation  temperature  is  higher  for  lower  yield,  T^  ~  y-0.032 
and  for  higher  altitude,  T^  ~  p~^^. 

For  sea  level,  for  Y  =  1,  and  for  p  =  5  bars  (cf.  Sec.  5e  for  this 
choice)  Brode's  calculations  give  T^  «  3*8;  then  (5-39)  yields  T^  =  1.08, 
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or  a  radiating  temperature  of  10,800°.  This  is  a  reasonable  result, 
though  appreciably  higher  than  the  temperatures  usually  observed.  How¬ 
ever,  as  we  shall  show  in  Secs.  :je  and  f,  our  calculation  gives  the  max¬ 
imum  temperature  reached,  and  is  likely  to  be  somewhat  too  high  due  to 
our  approximations. 

e.  Beginning  of  Strong  Cooling  Wave 

p 

Equation  (5*19)  gives  the  inward  speed  of  the  cooling  wave  in  gm/cni 
sec.  Precisely,  this  is  the  speed  at  which  the  point  of  temperature  T^ 
(radiating  temperature)  moves  relative  to  the  material,  once  the  cooling 
wave  is  fully  established.  But  even  if  there  is  no  cooling  wave ,  i.e., 
if  we  have  simply  adiabatic  expansion  behind  the  shock,  a  point  of  given 
temperature  T^  will  move  inward.  This  "adiabatic  motion"  is  the  minimum 
velocity  which  the  point  T^  can  have.  Therefore,  if  the  adiabatic  speed 
is  greater  than  (5*19)>  it  will  be  the  correct  velocity.  Of  course,  there 
will  still  be  a  cooling  wave  because  this  is  needed  to  supply  the  energy 
for  the  radiation;  this  "weak  cooling  wave"  will  be  described  in  Sec.  5g. 
But  its  inward  motion,  more  accurately  the  velocity  of  its  foot  (point  C 
in  Fig.  2),  will  not  be  determined  by  the  requirement  of  sufficient  energy 
flow,  (5.19) >  but  by  the  "adiabatic  speed"  which  we  shall  derive  from 
(3*15)*  Region  IV  of  Fig.  2  will  now  be  absent.  Thus,  outside  the 
cooling  wave,  at  point  C,  region  V  will  begin  immediately,  with  the 
temperature  distribution  given  by  (3*15) • 

It  is  therefore  important  to  determine  the  time  t  (and  pressure  p  ) 

El  3, 
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at  which  u  of  (5*19)  becomes  larger  than  the  adiabatic  speed  of  temper¬ 
ature  T_ .  Before  this  time  t  ,  we  have  a  weak  cooling  wave.  We  call 

X  cl 

this  Stage  C  I;  afterwards,  the  cooling  wave  is  strong,  and  the  radiation 
is  essentially  described  by  the  theory  of  Sec.  5d  (Stage  C  II).  To 
determine  the  point  of  separation  p  between  these  two  stages  is  a  refine- 

3 . 

ment  of  the  considerations  of  Sec.  Jd. 

The  temperature  in  the  adiabatically  expanding  material  behind  the 
shock  is  a  function  of  time  and  position,  given  by  hydrodynamics  and 
equation  of  state.  The  inward  motion  of  a  point  of  given  temperature 
relative  to  the  material  is  given  by 


-  P 


dr 

dt 


(9  log  T/St) 
0  log  T/p  3R), 


(5-40) 


where  the  subscript  r  means  that  the  partial  derivative  must  be  taken  at 
given  material  point  r,  not  at  given  geometrical  radius  R.  We  have  from 
(5-15) 


5  log  T  =  1.8 

p  SR  (7*  -  lJpR 


Similarly,  from  (5.28) 


o  5  log  p  _  1.2  3 

p  St  “  “  t 


0.8  7*  -  0.9 


(5-41) 


(5.42) 
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assuming  the  strong  shock  relation  p  ~  t"1  .  The  ratio  (5.40)  is  then 


D  ill  -  0 -9)(r*  -  l)  pH 

p  dt  175  71  ~ 


(5.43) 


Here  we  use  the  relation  (3.10) 


7*  -  1  H 


(5.44) 


and  find 


p£-£<7’  -0.9>p|  £ 


(5.45) 


where  H  has  been  labeled  ^  because  it  refers  to  the  radiating  tempera¬ 
ture.  Here  we  may  insert  (3.30)  and  (3.31),  viz.. 


5.  =  lj2  fl 

t  Z^Ap. 


(5-46) 


We  may  then  compare  the  result  with  (5.18),  the  velocity  of  the  cooling 


wave  (setting  JQ  =  0) 


Ho  -Hi 


(5.47) 


The  comparison  gives 


10n/2  ,  , 

9  K7 


3/2  H 

°'9)  Jp =  Jl 


(5.48) 
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2 

Inserting  (3.33)  for  J  ,  and  changing  the  unit  of  p  from  dynes/cm  to 
6  2 

bars  =  10  dynes /cm  we  get 


T ' 17/ 3 

9  1  j _ 

10^2  7-  0.9 


With  p0 


1.29  x  10  "5,  7' 


=  1.15  this  becomes 


P 


1 


(5.^9) 


Setting  now  p  =  p  and  choosing,  as  in  Sec.  5d,  T*  =  3*6  and  T  =  1.08 

X  U  C  X 

yields 

=  5*0  bars  (5.50) 

Thus  the  critical  pressure  is  5  bars,  which  was  the  reason  for  the  choice 

of  this  number  at  the  end  of  Sec.  5d.  There  it  was  shown  that  5  bars  and 

T1  =  3*6  leads  to  T,*  =  1.08  so  that  our  numbers  are  consistent, 
c  1 

We  had  to  rely  on  Brode’s  solution  to  find  T^  for  a  given  p;  this 
could  only  be  avoided  by  obtaining  an  analytic  solution  for  the  isothermal 
sphere  which  will  be  discussed  in  a  subsequent  paper.  Apart  from  this 
our  treatment  is  analytical,  making  use  of  the  equation  of  state  and  the 
absorption  characteristics  of  air. 
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f.  Maximum  Emission 


In  this  section  we  consider  Stage  C  I,  i.e.,  the  condition  when  the 
pressure  is  larger  than  p  ,  (5»50).  Then  the  radiating  surface  is  in 
the  adiabatic  region  described  in  Sec.  3c.  The  condition  for  its  posi¬ 
tion  is,  (3.21),  t  =  R/50.  From  Brode’s  curves,  the  position  of  a  point 
of  temperature  T  *  near  1  is  given  approximately  by 

R  =  0.78  Y1/5  t'-0*10  p-1/4  (5.51) 

This  expression  comes  purely  from  the  numerical  calculation,  except  that 
the  correct  dependence  on  yield  is  inserted.  Y  is  in  megatons,  p  in  bars, 
R  in  kilometers.  Equation  (5*5l)  holds  from  p  =  5  to  100  bars  within 
about  5$.  The  absorption  coefficient  is  given  in  (4.10),  which  yields 

1.1.  1  * 10-2  (M  T'-7 

“  po  W 


=  8  x  105  p"1*45 


cm 


(5.52) 


using  the  equation  of  state  (4.11).  Equating  this  to  1/50  of  (5.51) 
gives 


t’4.8 


4.5 


y-l/3  -1.20 


(5.53) 


T*  =  1.37  Y-0*07  p"0*25 


(5.54) 


Thus  the  radiating  temperature  increases  as  the  pressure  decreases. 


Since  p 


> 


s.  -1*2 
O 

T*~t°*5  (5.55) 

i.e.,  the  increase  of  temperature  with  time  is  fairly  fast.  This  is 
entirely  due  to  the  "opening  up"  of  the  shock,  i.e.,  the  decrease  in 
absorption  with  decreasing  density  (pressure) . 

For  given  pressure,  the  radiating  temperature  (5*5*0  is  slightly 
higher  for  smaller  yield,  an  effect  which  has  been  observed.  A  factor 
of  1000  in  yield  corresponds  to  a  factor  1.6  in  temperature,  hence  a 
factor  7  in  radiation  per  unit  area.  The  total  radiated  power  (assuming 
black  body)  at  given  p  is  proportional  to 

R2  t.U  „  *2/3-0.26  p-2/3-1  „  *0.40  p-5/3  (5.56) 

The  relatively  low  power  of  the  yield  is  remarkable  in  this  formula, 
which  describes  Stage  C  I.  Since  the  total  energy  radiated  should  be 
approximately  proportional  to  Y,  the  duration  of  Stage  C  I  is  then  pro¬ 
portional  to  Y^'^.  The  observed  time  to  the  second  radiation  maximum 

0.5  * 

is  about  proportional  to  Y  .  Our  theoretical  dependence  of  T  on  yield 

is  therefore  somewhat  too  strong.  The  time  dependence  of  (5*5 6)  is  quite 
strong,  about  as  t^. 

As  we  have  shown  in  Sec.  5e,  Stage  C  I  ends  when  the  pressure  reaches 
5  bars.  For  this  value  of  p,  and  for  Y  =  1,  (5*5*0  gives  T  =  0.92. 

This  is  slightly  less  than  the  T*  =  1.08  deduced  in  Sec.  5d  for  the  same 
p  and  Y  from  the  cooling  wave.  The  discrepancy  must  be  due  to  a  small 
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inconsistency  in  our  approximations. 


In  our  one -dimensional  theory,  the  end  of  Stage  C  I  marks  the  max¬ 
imum  of  radiation,  both  in  temperature  and  total  emission.  In  Stage  C  I 
the  radiating  temperature  increases  with  decreasing  pressure,  (5.54), 
because  the  material  becomes  more  transparent.  In  Stage  C  II  the  reverse 
is  the  case,  (5*39),  because  the  material  which  has  gone  through  the 
cooling  wave  becomes  thicker  with  time,  (5*32).  This  material  provides 
opacity  for  the  visible  light  from  the  fireball;  since  it  becomes  more 
opaque,  the  radiation  must  now  come  from  a  layer  of  smaller  absorption 
coefficient,  ( 5 * 35 ) >  and  therefore  of  lower  temperature.  As  (5*36)  shows, 
the  increase  of  thickness  (t  in  the  denominator)  is  more  important  than 

t/tc 

the  continued  decrease  of  density  (factor  p  '  J) .  These  results  will  be 
modified  in  the  three-dimensional  theory.  Sec.  6. 

The  maximum  temperature  has  been  calculated  as  T '  =  0.92  or  1.08, 
from  our  two  calculations;  it  is  clearly  close  to  1,  i.e.,  10,000°. 

This  number  is  not  too  much  out  of  line  with  observation  considering 
that  we  have  calculated  a  maximum.  In  fact,  the  transition  from  Stage 
C  I  to  C  II  cannot  be  sudden  as  we  have  assumed;  the  cooling  wave  must 
begin  gradually,  and  therefore  the  temperature  peak;  which  we  have  calcu¬ 
lated  will  actually  be  cut  off  (Fig.  3).  The  observable  maximum  may 
easily  be  1000°  lower  than  our  calculation. 
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Fig.  3.  The  calculated  increase  in  temperature  in  Stage  C  I, 
and  decrease  in  Stage  C  II  (straight  lines),  and  the 
estimated  actual  behavior. 


The  radius  of  the  radiating  surface  increases  with  time  in  Stage 
C  I.  Inserting  (5«510  into  (5*5l)  gives 


0.34  -0.225  .0.27 

R~Y"p  -'~t 


(5-57) 


In  Stage  C  II,  the  surface  moves  rather  rapidly  inward  relative  to  the 
material,  due  to  the  cooling  wave.  In  addition,  for  sea  level  explosions 
at  least,  the  pressure  is  no  longer  much  above  ambient,  so  that  the  out¬ 
ward  motion  of  the  material  slows  down.  Thus  the  geometric  radius  of  the 
radiating  surface  no  longer  increases  much,  and  soon  begins  to  decrease. 
Therefore,  the  total  radiation  will  reach  its  maximum  at  the  same  time 
as,  or  very  soon  after,  the  maximum  of  the  temperature. 
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g.  Weak  Cooling  Wave 


In  Stage  C  II  the  cooling  wave  dominates  the  radiation  (Sec.  5d) . 

In  Stage  C  I  the  cooling  wave  also  exists,  because  the  energy  of  the 
radiation  must  be  provided.  However,  the  speed  at  which  the  wave  proceeds 
inward  is  now  governed  by  the  adiabatic  expansion,  (5*^5) •  In  order  to 
obtain  the  correct  flux  of  radiation  at  the  radiating  surface,  we 
must  therefore  use  (5 .18)  in  reverse:  The  temperature  T^  at  the  inner 
edge  of  the  cooling  wave  (point  B  in  Fig.  2)  will  regulate  itself  in  such 
a  way  that  (5«l8)  is  satisfied,  with  u  given  by  (5«^+5)*  Using  (5 *46), 
we  thus  get  the  condition 


10  JK  (y' 

9 


H 


H0  -«! 


(J1  “  V 


(5-58) 


or  solving  for  and  inserting  numbers: 


HP.  I  ^  _ 1  =  _rP!_  JU  p-5/2 

H1  1  -  VJ1  1 


(5.59) 


Neglecting  J^,  we  find  that  H0  increases  rapidly  as  p  decreases. 

We  may  insert  (5*5J0;  then  the  right  hand  side  varies  as  p“-^2  ~  t^. 
Thus  the  cooling  wave  starts  very  weak  and  then  rapidly  increases  in 
strength.  Its  "head"  (point  B  in  Fig.  2)  is  at  first  close  to  its  foot 
(point  C).  As  time  goes  on,  it  moves  more  deeply  into  the  hot  material, 
eating  up  region  II  of  Fig.  2.  The  energy  which  is  made  available  for 
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radiation  is  essentially  the  difference  -  H^:  The  material  drops 
suddenly  in  temperature  from  TQ  to  as  the  cooling  wave  sweeps  over 
it,  and  the  energy  difference  is  set  free  for  radiation. 

The  term  depends  on  H^,  on  p,  and  on  the  temperature  gradient 
in  the  region  inside  the  cooling  wave  (region  II ) .  If  this  region  is 
adiabatic,  the  temperature  gradient  can  be  calculated  from  (3.15).  Since 
is  also  related  to  the  temperature  gradient  in  the  adiabatic  region 
IV,  the  ratio  Jq/J^  tends  to  unity  as  -*  H^.  Thus  the  left  hand  side 
of  (5*59)  will  have  a  certain  minimum  value.  This  seems  to  indicate 
that  there  is  no  cooling  wave  at  all  until  the  pressure  has  fallen  below 
a  certain  critical  value.  We  have  not  investigated  this  point  in  detail. 
It  is  possible  that  it  is  simply  related  to  the  break-away  of  the  lumi¬ 
nous  front  from  the  shock  wave,  i.e.,  the  beginning  of  Stage  C. 

Equation  (5*59)  describes  the  weak  cooling  wave  in  Stage  C  I  no 
matter  what  the  distribution  of  temperature  in  region  II.  The  stage 
comes  to  an  end  when  HQ  reaches  the  maximum  possible  value,  Hc«  There¬ 
after,  the  wave  described  by  (5*59)  becomes  inadequate  to  supply  the 
radiation  energy:  Since  the  enthalpy  can  no  longer  increase,  the  speed 
of  the  cooling  wave  has  to  increase.  This  speed  is  then  given  by  (5.19), 
and  Stage  C  II  has  begun. 
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6.  EFFECT  OF  THREE  DIMENSIONS 


a.  Initial  Conditions  and  Assumptions 

The  fact  that  the  fireball  is  spherical  causes  some  deviations  from 
the  one -dimensional  theory  of  Sec.  5*  In  this  section  we  shall  discuss 
Stage  C  II,  the  penetration  of  the  cooling  wave  into  the  isothermal 
sphere,  in  three  dimensions.  Initially,  we  assume  that  the  cooling  wave 
has  just  reached  the  isothermal  sphere;  we  call  the  corresponding  time 
t  =  t&,  and  the  pressure  is  the  critical  pressure  p^  =  5  bars  derived  in 
Sec.  5d»  Subsequently,  the  mass  of  the  isothermal  sphere  decreases  due 
to  the  cooling  wave. 

We  assume  that  the  material,  which  is  at  temperature  T  >  T^  =  4000° 

at  the  initial  time  t  ,  will  stay  in  this  temperature  range  throughout 

o. 

Stage  C  II.  This  is  reasonable  because  this  material  will  be  heated  by 
the  ultraviolet  radiation  coming  from  the  inside,  because  of  the  strong 
absorption  of  air  of  medium  temperature  (5000  to  6000°)  for  UV  (Sec.  4c). 
The  UV  heating  is  expected  to  compensate  approximately  the  cooling  due 
to  adiabatic  expansion  of  this  material;  this  is  confirmed  by  rough 
estimates  of  the  heating  and  cooling.  We  do  not  know  how  the  temperature 
is  distributed  in  the  ,(warm  layer"  between  the  radiating  temperature  T^ 
and  the  temperature  T^  =  4000°;  this  could  only  be  determined  by  a  detailed 
calculation  of  the  UV  radiation  flow  in  this  region.  We  assume  that  the 
distribution  is  smooth,  as  it  is  in  Stage  C  I,  and  that  therefore  the 
thickness  of  the  "warm  layer"  in  gm/cm^  is  directly  proportional  to  the 
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required  mean  free  path  of  visible  light  at  the  radiating  temperature 
2 

(also  in  gm/cm  )  which  in  turn  determines  itself. 

0 

We  take  the  initial  conditions  from  Brode,  using  his  curves  at  a 
time  when  the  inside  pressure  is  4.1  bars,  this  being  closest  to  p  =  5 
bars  of  al 1  the  curves  he  has  published.  This  corresponds  to  a  scaled 
(l-megaton)  time  t  =  0.325  sec.  At  this  time,  important  physical  quan¬ 
tities  are  as  given  in  Table  VIII  (dimensions  scaled  to  1  megaton).  The 
last  column  of  the  table  is  not  given  by  Brode  but  will  be  explained  in 
Sec.  6c.  In  the  table,  we  have  defined  a  quantity  proportional  to  the 
mass, 

m  =  p q1  J  p(R)R2  dR  (6.1) 

0 

where  R  is  measured  in  hundreds  of  meters.  It  is  interesting  that  the  mass 
of  the  warm  layer  is  much  (6.6  times)  larger  than  that  of  the  isothermal 
sphere.  Its  volume  is  about  2.4  times  larger  in  Brode’s  calculations. 

Table  VIII.  Conditions  When  Cooling  Wave  Reaches  Isothermal 
Sphere,  According  to  Brode 


Isothermal 

Warm  Layer 

Sphere 

Brode 

Sec.  6c 

Outer  radius,  meters 

385 

577 

438 

Average  density/ p^ 

0.79  x  10"2 

2.20  X  10-2 

11  x  10 

"Mass"  m  defined  in  (6.1) 

0.150 

0.99 

0.99 

Mean  temperature  T* 

3-0 

-- 

0.70 
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b.  Shrinkage  of  Isothermal  Sphere 

We  denote  the  "mass"  of  the  isothermal  sphere,  as  defined  by  (6.l), 
by  m^.  Then  this  mass  will  decrease,  due  to  the  progress  of  the  cooling 
wave  inward,  according  to 


JLR2 
dt  --p/l 


(6.2) 


Note  the  pQ  in  the  denominator  and  the  absence  of  the  factor  4- it,  both 

due  to  the  definition  (6.1).  The  density  of  the  isothermal  sphere  is 

nearly  uniform  and  will  be  denoted  by  p.  .  Initial  values  at  time  t 

is  a 

will  be  denoted  by  a  subscript  a.  The  isothermal  sphere  expands  adia- 
batically,  hence 


P 


is 


(6.3) 


and  therefore  at  any  time  t 


(6.4) 


The  speed  of  the  cooling  wave  is  given  by  (5*10),  thus 


u  = 


(6.5) 


where  JQ  has  been  set  equal  to  zero  because  there  is  no  appreciable  flow 


72 


of  radiation  inside  the  sphere.  is  the  enthalpy  inside  the  sphere, 
■which  decreases  adiabatically 


H. 


0 


(7-l)/7 


(6.6) 


is  the  enthalpy  at  the  radiating  surface.  It  will  be  shown  in  Sec.  6c 
that  the  temperature  is  nearly  constant  with  time,  so  that  in  good 

approximation  IL^  and  in  (6. 5)  are  constant.  We  shall  also  make  the 

poorer  approximation  that  «  H^.  Then  (6.5)  and  (6.6)  give 


u  =  u 

a 


(7-l)/7 


(6.7) 


Inserting  (6.4)  and  (6.7)  into  (6.2)  and  making  the  equation  dimen¬ 
sionless, 


(6.8) 


We  define 
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A 


u  t 
a  a 

P  R 
a  a 


(6.9) 


and  obtain  the  simple  differential  equation 


§y  _ 

dx 


Af 


(6.10) 


which  yields,  together  with  the  boundary  conditions, 


y  =  1  -  A 


X 

I 


f(x.)l-V37 


dx 


(6.11) 


Trie  constant  A  can  t>6  doteirnincd  from  See.  5^*  Tne  condition  tiiene 


is 


u 


a 


(6.12) 


where  p  is  the  density  outside  the  cooling  wave. 


According  to  (5*44) 


(6.13) 


Then  using  (f3 -43) 


u 

a 


ja  4  (y*  -  0.9)(7*  -  1)  Ra 
H  9  7T  t~ 

x  a 


(6.14) 
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Using  =  5 .6,  =  0.98,  H  ~  T5^5,  we  find  =  8.7.  Using  further 

7*  =  1.18,  A  in  (6.9)  becomes 

A  =  8.7  X  |  X  °*1'81^iB'<2,8  =  0,166  (6.15) 

which  is  a  pure  number,  and  small:  The  calculation  could  be  improved  by 
not  neglecting  compared  with  H^. 

According  to  (6.1l)  the  fireball  is  used  up  at  the  time  x  defined 
by 


dx' 


1 


(6.16) 


This  relation  is  valid  (granted  the  approximations  we  have  made)  what¬ 
ever  the  relation  between  pressure  and  time,  f(x).  This  generality  is 

useful  for  sea  level  explosions,  where  p  is  only  5  times  ambient  pressure: 

o. 

The  relation  between  p  and  t  can  then  be  taken  from  a  machine  calculation 
(or  ob s ervat ion ) . 

For  higher  altitude,  let  us  say  h  >  10  km,  p  =  5(p,/pr,)^^  is 
sufficiently  above  the  ambient  pressure,  P^/p q>  so  that  the  shock  is 
still  strong  and  we  may  use 

f (x)  =  x1*2  (6.17) 

Then  (6.1l)  becomes 
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(6.18) 


y  =  1  -  A 


x2.2-0.k/y  _  ± 

2.2  -  0.4/7 


From  7  =  1.18, 


y  = 


12.2 


1.86 

X 


11.2 


(6.19) 


and  the  isothermal  sphere  disappears  for 


x2  =  12. 21/1*86  =  3.84 


(6.20) 


For  these  higher  altitudes,  then,  the  time  when  the  isothermal  sphere 
disappears  is  a  fixed  multiple  of  the  time  when  it  is  first  reached  by 
the  cooling  wave.  This  multiple  depends  only  on  y}  and  on  the  ratio 
\/\  internal  to  external  enthalpy  at  time  t&.  The  pressure  at  the 
time  x2  is 


—  fef 


(6.21) 


For  sea  level,  f(x)  increases  more  slowly  (the  pressure  decreases 
more  slowly)  with  time;  hence  it  takes  somewhat  longer  to  use  up  the  iso¬ 
thermal  sphere.  Conversely,  the  pressure  at  the  time  t„  =  t  x  will  have 

2  a  a 

decreased  by  a  smaller  factor  from  p  . 

a 

For  the  simple  case  of  higher  altitude,  we  can  use  (6.19)  to  calcu- 
late  the  fraction  of  the  mass  Y  which  will  still  be  in  the  isothermal 
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sphere,  its  radius  R^/R^,  and  other  physical  data.  Some  of  these  are 
given  in  Table  IX.  It  is  seen  from  the  table  that  the  mass  decreases, 
first  fairly  uniformly  and  rapidly  (as  if  it  would  go  to  zero  at  x  =  2.8), 
then  more  slowly  (because  it  is  proportional  to  y^),  while  the  radius 
first  expands  slightly,  then  shrinks  slowly  and  at  the  end  very  rapidly. 
The  latter  phenomenon,  the  rapid  shrinking  of  the  apparent  fireball,  may 
not  be  observable  because  the  bomb  debris  becomes  visible  and  is  likely 
to  be  still  opaque.  But  some  shrinkage  of  the  fireball  should  be  open 
to  observation. 


Table  IX.  Development  of  Isothermal  Sphere  and  Warm  Layer 
due  to  Cooling  Wave 


X  =  t/ta 

1.0 

1.2 

1.5 

2.0 

2.5 

3.0 

3.5 

5.84 

1 

0.964 

0.899 

0.766 

0.599 

0.402 

0.188 

0 

II 

1? 

1 

0.895 

0.726 

0.450 

0.215 

O.O65 

0.0066 

0 

VRa 

1 

1.029 

1.051 

0.970 

0.818 

0.584 

0.287 

0 

VRa 

1.137 

1.194 

1.248 

1.296 

1.526 

1.562 

1.426 

1.485 

L* 

0.137 

0.129 

0.126 

0.129 

0.150 

0.180 

0.214 

0.247 

The  last  two  lines  of  Table  IX  will  be  explained  in  Sec.  6c. 
c.  The  Warm  Layer 

We  want  to  assume  that  the  outer  edge  of  the  warm  layer,  the  4000° 
temperature  level,  stays  fixed  in  the  material  (Sec.  6a).  We  wish  to 
calculate  the  thickness  of  the  warm  layer  in  gm/cm^. 
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p(R)  cLR 


(6.22) 


where  Rg  is  the  outer  edge  of  the  layer.  The  original  mass  of  the  warn 
layer  m^  is  much  larger  than  that  of  the  isothennal  sphere 

VKi  =  B  =  6,6  (6.23) 

At  a  later  time,  the  mass  of  the  warm  layer  is  then 

“2  =  mjB  +  1  -  y5)  (6.24) 

* 

It  does  not  change  much. 

The  temperature  in  the  warm  layer  goes  from  10,000  to  4000°.  For 
simplicity  we  assume  that  the  density  corresponds  to  the  average  temper¬ 
ature  of  7000°.  Then,  at  the  initial  pressure  used  in  Table  VIII,  p  = 

cl 

4.1  bars,  the  density  of  the  warm  layer  is  p  =  11  X  io“2  p  =  l4p  . 

v  0  a 

This  is  much  higher  than  in  Brqde's  calculations,  Table  VIII:  In  his 
calculations,  the  cooling  wave  has  not  reached  the  isothermal  sphere; 
with  our  assumptions  it  has.  The  difference  is  due  to  the  different 
opacities  assumed;  it  has  the  consequence  that  our  warm  layer  is  geo¬ 
metrically  much  thinner  than  his  (Table  VIII,  last  column).  Subsequently 
the  density  decreases  with  pressure,  not  adiabatically  but  isothermally; 
in  accord  with  Gilmore's  formula  (3.13),  we  assume 
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(6.25) 


with  C  =  14. 

On  this  basis  we  calculate  the  outer  radius  Rg  of  the  warm  layer 
and  find 

R^  =  R:?  +  R5  B  +  1h~  ^  Cf(x)]10^  (6.26) 

d  JL  cl 


Using  (6.17)  we  can  then  calculate  R„/R  .  We  give  this  quantity  in  the 
second  last  line  of  Table  IX.  In  the  last  line,  we  have  given  the  mate¬ 
rial  thickness  of  the  warm  layer,  in  relative  units. 


L 


I 


(6.27) 


It  is  seen  from  the  table  that  Rg  increases  only  slowly;  the  main 
change  in  Rg  -  R^  at  later  times  is  therefore  due  to  the  decrease  of  R^. 
The  material  thickness  l'  first  decreases  very  slightly;  this  continues 
about  as  long  as  R^  increases,  and  is  due  to  the  fact  that  about  the  same 
mass  of  warm  material  gets  distributed  over  a  larger  area.  Later  on,  L* 
increases  while  R^  decreases.  Until  x  =  2.5  the  change  of  L*  remains 
less  than  10$,  and  after  x  =  2.5  the  calculation  is  probably  meaningless 
because  the  bomb  debris  comes  into  view.  Therefore,  we  may  assume  L' 
constant,  and  thus  the  optical  mass  absorption  coefficient  (in  cm  /gm) 
at  the  radiating  layer  will  also  be  constant.  In  contrast  with  this 
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result  of  the  three-dimensional  considerations,  the  one -dimensional  for¬ 
mula  (5*35)  yields  a  decrease  of  p/p  as  t~^. 

Using  (5*33) >  a  constant  mass  absorption  coefficient  requires 


~  p_1/5n  ~  p“1//20 


(6.28) 


Thus  in  three  dimensions,  the  temperature  continues  to  increase  slowly 
after  the  cooling  wave  reaches  the  isothermal  sphere.  If  we  assume  that 
times  up  to  x  =  2.5  in  Table  IX  are  significant,  the  pressure  decreases 
by  about  a  factor  of  3  in  Stage  C  II,  and  the  temperature  of  the  radiating 
layer  increases  by  about  5$  according  to  (6.28).  The  radiated  power, 
being  proportional  to  R^  T^  ,  may  increase  by  about  10$  up  to  x  =  2, 
and  then  decreases  due  to  the  shrinkage  of  the  radiating  surface. 

One  prediction  of  this  theory  is  that  the  temperature,  as  well  as 
the  radiated  power,  has  a  rather  flat  second  maximum  while  the  isothermal 
sphere  radiates  away  its  energy.  The  slow  variation  of  T  justifies  the 
treatment  of  J  as  constant  in  Sec.  6b. 

7.  TRANSPARENT  FIREBALL  (Stage  D) 

After  the  isotnermal  spnere  has  been  eliminated  by  the  cooling  wave, 
the  fireball  is  transparent  if  we  neglect  the  effect  of  the  bomb  debris. 

At  this  time  the  pressure  is  given  by  (6.21),  the  average  temperature  of 
the  Varm  region  is  about  7000°,  and  the  corresponding  density  is  about 
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(7.1) 


p 

—  =  2.3  X  10 
P0 


0.37 


If  we  assume  that  the  radius  R  of  the  isothermal  sphere  scales  with 

cl 

-I/? 

J  (which  may  be  wrong)  the  radius  of  the  warm  sphere  is  now 


Rb  =  570 


(7-2) 


where  Y  is  in  megatons.  Therefore,  along  a  radius,  the  amount  of  warm 
material  is  about 


PyRjj  =  1*7  Y1//5  gm/cm2  (7.3) 

almost  independent  of  ambient  air  density. 

Taking  P^/Pq  =  10 "2  (which  corresponds  to  P^/Pq  10”^)  and  an 
average  temperature  of  7000°,  Meyerott  *s  tables  give  for  the  visible: 

cm2/gm  (7.4) 

so  that  the  radius  represents  0.7  optical  mean  free  path  at  1  megaton. 
This  is  nearly  transparent,  and  the  fireball  becomes  rapidly  more  trans¬ 
parent  as  it  cools  down  by  further  emission  of  radiation.  The  emission 
of  radiation  is  then  proportional  to  the  opacity;  for  each  material 
element. 
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(7-5) 


3?-  «=  -  teaT4 
dt 

p 

where  a  is  the  Stefan-Boltzmann  constant  and  K  the  opacity  in  cm/gm, 
i.e.,  the  absorption  coefficient  averaged  over  a  Planck  spectrum  (dis¬ 
regarding  spectral  regions  which  are  still  black) .  Since  K  increases 
rapidly  with  temperature,  the  hottest  region  near  the  center  will  cool 
fastest,  so  that  the  temperature  tends  to  become  more  uniform. 

For  the  same  reason,  the  radiative  cooling  will  effectively  stop  at 

a  temperature  of  about  5000°.  At  this  temperature  and  p  =  O.Olp^,  it 

takes  about  5  seconds  to  cool  the  fireball  by  10$.  At  4000°,  this  takes 

25 

about  200  seconds.  Gilmore  Gias  calculated  curves  of  cooling  times  for 
transparent  bodies  at  various  densities  as  a  function  of  the  final  temper 
ature. 

Depending  on  the  ambient  density,  the  pressure  may  or  may  not  have 
decreased  to  ambient  pressure  when  the  temperature  has  decreased  to  5000° 
Even  if  it  has  not,  the  subsequent  adiabatic  expansion  will  not  lower  the 
temperature  much  further.  (Our  theory  is  not  applicable  to  very  low 
ambient  densities  because  there  the  isolation  of  isothermal  sphere  from 
the  outside  never  takes  place.  Probably  the  limit  of  applicability  is 
about  P^/pQ  =  10  Therefore,  and  because  of  the  small  value  of  y  -  1, 
expansion  cannot  be  large.)  Therefore,  after  both  radiation  and  hydro¬ 
dynamics  have  effectively  stopped,  the  fireball  is  left  at  a  temperature 
not  much  below  5000°. 

Any  further  cooling  can  only  be  achieved  by  the  rise  of  the  fireball 
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due  to  its  bouyancy,  and  the  turbulent  mixing  associated  with  this  rise. 

This  is  a  slow  process,  taking  tens  of  seconds. 

Since  the  emission  is  now  proportional  to  the  absorption  coefficient, 

the  molecular  bands  will  now  appear  in  emission  while  in  earlier  stages 

they  appear  in  absorption.  This  has  been  observed. 

The  debris,  at  the  center  of  the  fireball,  contains  metals  and 

therefore  is  likely  to  be  opaque  at  lower  temperatures.  Therefore  the 

debris  may  well  be  opaque  after  all  the  air  has  become  transparent.  The 

debris  usually  has  a  ragged  shape  due  to  Taylor  instability.  Recently, 

26 

Longmire  has  given  a  tentative,  quantitative  theory  of  this  instability 
in  debris  expansion.  Because  of  its  higher  opacity,  the  debris  may  cool 
to  a  lower  temperature  than  the  surrounding  air. 
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